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SYMBOLS  AND  TERMINOLOGY 


c,m  = subscripts  denoting  'carrier'  and  'modulation' 
f = frequency 

g = gravitational  acceleration 

= empirical  cosine  coefficients  of  probe  angular  response 
k = 2n/L  = wavenumber 
n = probe  response  exponent 

n^  = radian  frequency  of  biharmonic  wave  components 
t = time  referred  to  arbitrary  origin 
u = horizontal  component  of  u 
u = fluid  flow  velocity 
w = channel  width 

x = coordinate  distance  from  the  paddle  mid-position  positive  along  the 
channel 

y * cross-channel  coordinate 

z * vertical  coordinate  referred  to  still  water  level. 

a = asymmetry  coefficient,  also  paddle  servo  voltage,  probe  directional 
coefficient 

B=(kn  - k r>  )/k  n = biharmonic  wave  parameter 
v m c c nr  c c r 

C = phase  velocity 

C = group  velocity 
S 

Cj  = ith  order  Stokes  phase  velocity 
C*  = Stokes  limiting  wave  phase  velocity 
D = channel  depth 
E = total  energy  density 
Ej.  = kinetic  energy  density 

Ep  = root-mean- square  surface  elevation  (potential  energy  density  = PgEp) 
Epx  = spatial  equivalent  of  Ep 
F = potential  energy  flux 
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F0  * Fx(G  = 0.19)  = normalizing  energy  flux 

Ft  = total  energy  flux 

F(t)  = initial  growth  function 

G = -Lji(3w/3x)/w  = growth  rate 

H =•  wave  height,  crest  to  trough 

H/L  = wave  steepness 

L = wavelength 

Lj  = q/f 

L*  = C*/f  = length  of  Stokes  limiting  wave 
Q = total  energy  loss  rate 
= breaking  energy  loss  rate 
Qr  = reflected  energy  loss  rate 
Qy  = viscous  energy  loss  rate 
T = 1/f  - wave  period 
U = Ouid  flow  speed 

V = anemometer  signal  voltage 

* 

•» 

VQ  " st:.:l  water  signal  voltage 
a = n^/n,  = biharmonic  frequency  ratio 

t l 

A = finite  time  or  space  increment 
n = surface  elevation  referred  to  still  water 
6 = angle  of  flow  with  respect  to  hot-film  probe  axis 
v = kinematic  viscosity  of  water 
p = water  density 

4>  = flow  angle  with  respect  to  orthogonal  probe  bisector 
u = 2nf  = wave  radian  frequency 

= overbar  indicating  ensemble  average 
* = subscript  for  Stokes  limiting  wave 


-*>  = vector  sign 

V = 3/3x  + 3/3y  * differential  operator 


ABSTRACT 

This  study  represents  the  first  known  attempt  to  examine  in  some 
detail  the  flow  dynamics  of  the  terminal  growth  and  breaking  of  periodic 
surface  waves  in  deep  water,  under  controlled  and  reproducible  laboratory 
conditions.  This  work  was  considered  a necessary  preliminary  to  future 
studies  of  storm  wave  breaking,  a phenomenon  of  considerable  engineering 
importance  that  has  received  little  scientific  attention. 

In  the  present  study,  a tape-controlled  servo-paddle  system  pro- 
duced single-frequency,  deep  watf/r  wave  trains,  which  were  laterally  con- 

tf 

verged  to  the  breaking  point  within  a tapered  channel.  Using  high-response 
elevation  and  velocity  sensors profiles  and  internal  velocity  fields  of 
growing  waves  were  determined  as  functions  of  time  and  distance  before 
and  after  breaking.  Parametric  analysis  of  these  data  yielded  the  following 
descriptions  of  consecutive  growth  stages. 


1)  For  wave  steepnesses  smaller  than  0.1,  wave  profiles  remained 
fairly  symmetric,  and  all  observed  prope»ies  corresponded 
reasonably  to  those  given  by  Stokes  fifth -order  theory. 

* 

2)  For  steepness  within  the  rai'ge  0.1<H/L<0.12,  wave  development 
was  characterized  by  increasing  profile  asymmetry,  steepening 

of  the  forward  face,  and  forward*  tilting  a'nd  upward  concentration 
of  internal  isolines  of  const* t .velocity.  These  changes  were 
confined  to  the  upper  30  percent  of  the  wave  height.  The  rate 
of  increase  of  potential  energy  density  along  the  channel 
before  breaking  was  found  to  be  inversely  related  to  channel 
width,  consistent  with  conservation  of  energy  flux. 

3)  Wave  breaking  occured  upon  attainment  of  steepnesses  within  the 
range  0. 125<H/L<0. 140,  and  was  associated  with  near-verticality 
of  the  forward  crest  face,  and  distortion  of  velocity  isolines 
so  as  to  wrap  around  a jet  issuing  from  the  forward  face  just 
below  the  crest.  The  vertical  profiles  of  maximum  horizontal 
velocity  were  essentially  hyperbolic,  and  maximum  jet  velocity 
exceeded  the  local  crest  velocity  by  as  much  as  10  percent. 

The  latter  velocity  was  found  to  closely  approach  the  theoretical 
limit,  C*=  1.2  g/2wf,  for  Stokes  waves  of  meximum  height. 

4)  Breaking  intensity  • as  determined  from  potential  energy  loss 
rate  - was  found  to  be  correlated  with  the  rate  of  potential 
energy  increase  before  breaking,  indicating  that  any  determin- 
istic theory  for  breaking  must  include  some  prescription  of  a 
wave’s  previous  history. 

5)  Post-breaking  energy  loss  rates  were  found  to  be  roughly  in 
equilibrium  with  energy  input  resulting  from  continuing  channel 
convergence. 
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1.  INTRODUCTION 


The  present  experiments  were  initially  proposed  as  the  first  phase 
of  a continuing  investigation  of  the  breaking  of  deep  water  waves,  such 
as  occurs  in  Nature  on  all  scales,  from  incipient  whitecapping  under 
light  winds  to  massive  plunging  of  higher  waves  under  storm  conditions. 
The  principal  argument  for  the  overall  program  is  that  wave  breaking  is 
very  little  understood,  but  of  manifest  theoretical  and  practical  impor- 
tance to  many  areas  of  ocean  science  and  technology. 


o It  is  generally  conceded  to  be  the  principal  energy  loss  mech.a".ism 
regulating  the  growth  of  sea  state  toward  spectral  equilibrium 
under  steady  winds  (Phillips,  1969,  p.  109), 

o It  is  conceptually  a more  efficient  mechanism  for  transferring 
momentum  from  wind  to  mean  surface  flow  than  viscous  dissipation 
by  non-breaking  waves  (Longuet-Higgins,  1969a),  or  direct  tang- 
ential stress  (Snyder  & Cox,  1966). 

o At  advanced  sea  states,  breaking  waves  exert  by  far  the  largest 
unitized  forces  against  ships  and  interfacial  structures,  with 
all  their  implications  to  ship  design  and  operation. 

In  view  of  the  above,  one  is  led  to  wonder  why  almost  all  wave- 
oriented  research  within  the  past  two  decades  has  been  directed  towards 
wave  growth  mechanisms,  3s  opposed  to  wave  breaking.  There  seem  to  be 
at  leust  two  reasors.  Wave  break/ng — a-.jdefined  by  turbulent  energy  loss- 
-is  a non-stationary  and  nonlinear  process  for  which  there  is  presently  no 
adequate  mathematical  description.  Second,  •->  sn  under  steady  winds,  the 
sea  surface  is  highly  irregular  and  non-repeatable  in  time  or  space,  so 
that  its  properties  aie  not  readily  definable  on  a wave-by-wave  basis. 
Instead,  it  has  been  found  more  practical  to  work  with  statistical  or 
spectral  representations  that  define  only  average  properties  and  their 
probabilistic  extremes. 

In  a macroscopic  sense,  however,  breaking  i_s  a wave-by-wave  pheno- 
menon, which  is  believed  to  occur  sporadically  whenever  random  spectral 
components  constructively  interfere  to  locally  increase  the  height  of  a 
particular  wave  to  the  point  where  it  becomes  unstable.  While  the  criteria 
for  wave  instability  are  not  known *,  where  the  surface  spectrum  is  reason- 
ably narrow,  as  manifested  by  fairly  regular  wave  groups,  Donelan,  et  al . 

^Phillips  (1969)  supposes,  by  analogy  to  the  well  known  stability  limits 
for  progressive  and  standing  gravity  waves,  that  the  probability  of 
irregular  waves  breaking  is  some  function  of  the  local  fluid  acceleration 
at  the  surface  as  a fraction  of  the  gravitational  acceleration,  g. 
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(1972)  have  observed  (and  confirmed  by  simple  experiment)  that  wave  breaking 
occurs  systematically  at  the  group  antinodes,  and  at  intervals  equal  to 
roughly  twice  the  period  of  the  principle  waves  present.  This  result, 
coupled  with  increasing  evidence  (Rye,  et  al . , 1974)  that  storm  spectra  are 
much  narrower  than  those  proposed  some  years  ago  by  Pierson  and  Moskowitz 
(1964),  suggest  the  possibility  that  a careful  series  of  laboratory  experi- 
ments might  yield  descriptions  of  the  breaking  process  that  have  applica- 
tion to  the  real  ocean. 

While  interactive  breaking  is  readily  produced  without  wind  in  the 
laboratory,  and  is  currently  proposed  for  the  next  phase  of  this  investi- 
gation, it  was  considered  desireable  to  first  explore  the  growth  toward 
breaking  of  single- frequency,  uniform  wave  trains  by  slow  lateral  conver- 
gence within  a tapered  channel  whose  depth  was  great  enough  not  to  signi- 
ficantly affect  the  motion.  This  method  of  approach  had  the  following 
advantages : 

o wave  growth  was  fully  controllable  and  reproducible,  and  permitted 
detailed  observation  of  wave  properties  at  all  stages, 

o it  most  nearly  meets  the  assumptions  of  steady-state  theories  for 
waves  of  finite  height, 

o it  permitted  control  of  breaking  intensity  by  varying  only  initial 
wave  amplitude,  although  frequency  was  also  varied  to  give  a 
greater  range  of  initial  wave  steepness. 

The  basic  research  objectives  of  this  study  include: 

1)  for  representative  initial  wave  steepnesses,  to  make  sufficiently 
detailed  measurements  of  surface  elevation  and  internal  fluid 
velocities,  as  functions  of  time  and  position,  so  as  to  reason- 
ably define  the  flow  fields  within  waves  approaching  the  breaking 
point, 

2)  to  compare  these  results  with  relevant  theories, 

3)  to  attempt  to  parameterize  breaking  and  breaking  intensity 
(energ/  loss)  in  terms  of  observable  variables,  as  a guide  for 
future  laboratory  or  field  experiments  with  interactive  breaking. 


2.  BACKGROUNP 


The  present  investig 'tion  was  closely  patterned  after  a similar  set- 
of  convergent  breaking  experiments  reported  by  Lc  Mchaute  ct  al . (1968). 
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The  latter  were  confined  to  shallow  water  waves,  were  conducted  in  a 
channel  about  half  the  present  dimensions,  and  were  limited  to  waves  about 
one- sixth  as  high  (10  cm).  Data  consisted  of  stroboscopic  photos  of  sur- 
face profiles  and  the  orbital  motions  of  neutral-density  particles  of 
nitrile  rubber.  T!‘3ir  methodology  and  results  are  of  interest  here  for 
several  reasons.  First,  the  experiment  was  designed  to  most  nearly  ideal- 
ize the  theoretical  assumption.;  of  steady  motions  in  uniform  depth. 

Second,  the  observed  distributions  c horizontal  particle  velocity  beneath 
the  wave  crests  (considered  *-o  b-  the  most  sensitive  test  of  theoretical 
validity)  were  inconsistent  w.t vnuse  predicted  by  any  of  eight  theories 
examined.  Tne1.  included  shallow-water  expansions  of  Stokes  waves  from 
first-  to  1 •' fth-erder,  C ...al  waves,  and  Solitary  waves.  Within  its 
region  of  5>  'r . inability , cue  Stokes  fifth  order  solution  gave  the  best  fit 
to  the  obse /ved  wave*  pit files,  and  vertical  velocity  distributions  of 
smaller  n n-b^eaxing  ’..ares,  although  the  observed  velocity  magnitudes  were 
smaller  l”  n 25  per  cr . ■ . T'v  silence  is  that  no  present  theory  reason- 
ably predicts  the  f.’ow  fields  -ithin  steep  shallow-water  waves. 

In  deep  water,  the  bottom  boundary  condition  for  irrotational  waves 
becomes  much  simpler,  and  o.ie  is  left,  effectively,  with  Stokes'  solutions 
of  higher  order  for  waves  of  fin.te  height  (Kinsman,  1965).  Partly 
because  of  their  complexity  and  partly  because  stochastic  sea  state  models 
are  consequently  li"-‘+cd  to  linear  combinations  of  small -amplitude  waves, 
these  solutions  have  not  received  much  practical  application.  The  current 
situation  is  summarized  by  Mr :.xmo>er  and  Kutzbach  (1965),  who  also  present 
a numerical  method  of  ootaining  solutions  to  any  order,  although  the  condi- 
tions of  profile  symmetry  and  steady  motion  are  still  retained.  As  later 
shown,  their  calculated  wave  profile  for  a 15th-order  wave  of  steepness, 
H/L  = 0.1,  is  in  good  agreement  with  those  observed  in  these  experiments 
during  the  earliest  growth  stage. 

The  limiting  form  of  a Stokes  wave  has  been  studied  extensively 
(Lamb,  1932,  pp.  418-419).  The  development  is  based  on  the  kinematic 
assumption  that  the  flow  velocity  at  the  crest  cannot  exceed  the  phase 
speed,  and  well  known  properties  include: 

o a cuspated  crest  symmetric  about  the  vertical,  with  an  included 
angle  of  120°; 

o a phase  velocity  20  percent  greater  than  that  for  infinitesimal 
waves  of  the  same  frequency:  C*  = 1.2g/27if; 
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o a maximum  height  equal  to  14.2  percent  of  the  maximum  wavelength, 

L*  = CJ  f; 

o a maximum  crest  acceleration  of  0.5g  directed  interiorly  (Longuet- 
Higgins,  1969b) . 

Various  recent  theoretical  attempts  to  describe  breaking  all  start  from  the 
Stokes  limiting  form  and  are  restricted  to  steady  motion  1 (Price,  1971; 
Longuet-Higgins,  1973,  1974).  In  contrast,  our  results  indicate  that 
breaking  is  strongly  influenced  by  the  rate  at  which  waves  grow,  so  that 
any  realistic  description  of  breaking  must  include  some  specification  of 
growth  rate. 


3.  FACILITIES  AND  INSTRUMENTATION 
3 . 1 Wave  Channel  and  Convergent  Barrier 

All  experiments  were  conducted  in  the  2.4x2.4x43m(8x8x  150 
ft)  wind-wave  channel  in  the  Scripps  Institution's  Hydraulic  Facility  (Fig. 
1).  For  the  convergent  breaking  experiments,  a rigid  barrier  2.4  m high 
and  24.3  m long  was  constructed  of  heavy  plywood  backed  by  steel  channel 
frames.  The  frames  were  bolted  together  and  through  the  channel  bottom, 
and  were  laterally  braced  at  the  top  and  mid-height  so  that  there  were 
negligible  deflections  during  wave  passage.  The  barrier  was  faired  into 
one  channel  wall  about  1 m from  the  wave  paddle  and  converged  toward  the 
opposite  wall  at  a rate  dw(x)/dx  ? 0.10  for  about  19  m,  where  it  terminated 
in  a hingepoint  at  a convergence  width,  w = 51.5  cm.  An  additional  6.1  m 
hinged  barrier  extension  was  designed  to  be  adjustable  over  the  range 
0 < dw/dx  < 0.05,  and  could  be  rigidly  secured  in  any  position. 

It  was  originally  hoped  that  the  hinged  barrier  might  be  adjusted  by 
trial  so  as  to  maintain  breakers  of  constant  height,  but  this  was  found 
impossible  because  of  strong  reflections  from  the  1:8  beach  slope  just 
beyond  the  convergence.  After  some  experimenting,  it  was  determined  that 
reflections  could  be  virtually  eliminated  by  closing  the  convergence  to  a 
terminal  gap  of  30  cm,  adding  a divergent  extension  at  an  angle  of  about 
45°  across  the  slope,  and  covering  the  upper  part  of  the  slope  with  wire 

^A  necessary  and  sufficient  condition  for  steady  motion  that  is  that  flow 
velocities  are  a function  of  phase  only,  <j>  = kx  - 2nft. 
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cages  of  metal  turnings.  All  reported  observations  were  ultimately  conducted 
within  the  1/10  convergent  section,  and  a 1.2  m plexiglass  window  was  cut 
into  the  channel  wall  18  m from  the  wave  paddle  for  sidewall  photography. 

Figure  1 also  shows  the  coordinate  system  used  throughout  these 
experiments;  the  origin  was  located  at  the  intersection  of  the  still  water 
level  and  the  paddle  mid-position,  with  the  x-axis  positive  away  from  the 
paddle,  and  the  z-axis  positive  upwards.  While  the  local  convergence  width, 
w(x),  remained  constant,  the  water  depth,  D,  was  varied  with  wave  frequency, 
as  later  described. 
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3 . 2 Wave  Generation 

The  wave  generator  is  an  integral  part  of  the  channel  complex  and  was 
used  without  modification.  It  consists  of  a 2.5  cm  plywood  paddle  rein- 
forced by  tapered  aluminum  backing  struts,  and  dependent  from  a double- 
articulated  trapezoidal  framework  bolted  to  the  building  roof  structure 
(Fig.  1).  By  changing  the  spacing  of  the  lower  hinge  points,  the  paddle 
motion  can  be  adjusted  to  pivot  about  any  horizontal  axis  from  infinitely 
far  below  the  channel  (planar  piston)  to  above  the  channel  top.  For  these 
experiments,  the  paddle  was  set  to  pivot  at  the  channel  bottom,  which  posi- 
tion was  found  to  produce  deep  water  waves  whose  spectra  were  essentially 
free  from  unexpected  energy  peaks  (App.  B) . To  minimize  cross-channel 
reflections  that  excited  side  modes,  the  paddle  was  fitted  with  a canted 
front  face  whose  normal  bisected  t lie  convergence  angle.  This  resulted  in 
very  uniform  wave  crests  which  sloped  up  barely  perceptibly  at  the  sidewalls 
(Fig.  5),  causing  the  initiation  of  breaking  at  these  points  an  instant 
before  it  occurred  along  the  remainder  of  the  crest. 

The  wave  paddle  was  driven  by  a hydraulic  piston  whose  position  was 
controlled  by  a potcntiometric  hydraulic  servo  system  slaved  to  an  electri- 
cal analog  input  voltage.  Maximum  paddle  excursions  corresponded  to  a 
voltage  swing  of  about  ±7  volts,  and  0-volts  represented  the  paddle  mid- 
position. The  actual  wavemaking  instructions  consisted  of  computer- 
produced  analog  signals  recorded  on  FM  tape.  These  were  demodulated  and 
played  back  to  the  servo  to  generate  completely  reproducible  wave  series  as 
often  as  desired,  the  wave  amplitude  being  controlled  by  the  demodulator 
setting. 
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Because  the  initial  wave  coming  off  the  paddle  traveled  at  roughly 
the  limiting  phase  speed  in  shallow  water,  whereas  the  energy  ’front* 
moved  essentially  at  deep  water  group  velocity,  it  was  necessary  to  smooth- 
ly increase  the  paddle  stroke  from  zero  amplitude  to  its  asymptotic  uniform 
value  in  order  to  prevent  premature  breaking  of  those  waves  preceding  the 
energy  front.  This  was  accomplished  by  incorporating  an  initial  growth 
function,  F(t),  in  the  paddle  displacement  equation 

A = AQF(t)  cos (2mft  + b) , (1) 

where  A and  Aq  are  the  input  and  maximum  servo  voltages,  respectively,  and 
b is  an  arbitrary  constant.  Typical  input  forms  for  the  function,  F(t)  = 
tanh(2irft/n  + 9/4n) , are  shown  in  Figure  2.  For  each  combination  of  f and 
nQ,  the  constant  n was  adjusted  by  trial  to  obtain  a continuous  succession 
of  uniform  breakers  at  a fixed  position  in  the  channel.  In  most  cases,  the 
wave  series  was  terminated  by  a ramp  function  after  about  30  seconds.  Scale 
and  time  reference  to  the  paddle  motion  were  produced  by  a motion  transducer 
mounted  on  the  paddle  arm. 


3.3  Velocity  Measurements 

Utilizing  an  orthogonal  pair  of  DISA  55A81  wedge-shaped  hot  film 
probes  coupled  to  DISA  55I|01  constant  temperature  anemometers,  over 
600,090  individual  determinations  of  flow  speed  and  direction  within  gro.- 
ing  waves  were  made  during  the  course  of  these  experiments.  The  probe 
shanks  were  encased  in  a watertight  aluminum  turret  with  their  emergent 
tips  abutting  at  right  angles  (Fig.  3).  Since  the  anemometer  output  from 
a single  probe  is  an  analog  voltage  that  is  functionally  related  to  the 
product  of  a velocity  term  and  a direction  term,  the  signals  from  two 
orthogonal  probes  whose  calibration  constants  are  known  provide  an  unambig- 
uous measurement  of  the  flow  vector  for  two-dimensional  flow  in  the  plane 
defined  by  the  probe  axes.  While  the  resolution  of  flow  speed  in  this  plane 
is  practically  limited  by  backflow  to  the  90°  angle  subtended  by  the  probe 
axes,  this  window  could  be  broadened  by  rotating  the  turret  in  45°  incre- 
ments in  either  direction. 
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The  probe  turret  was  mounted  on  a vertical  support  strut  dependent 
below  an  instrument  cart  that  could  be  rolled  to  any  position  along  the 
wave  channel  (Fig.  4).  The  turret  plane  was  parallel  to--and  25  cm 
distant  from — the  channel  wall,  or  about  in  the  middle  of  the  convergence 
at  its  hinge  point.  Because  principal  interest  was  directed  toward  meas- 
uring high  crest  velocities  having  components  in  the  direction  of  wave 
propagation,  and  because  negative  trough  velocities  were  generally  too  low 
to  resolve  accurately,  all  measurements  were  made  with  the  probes  oriented 
in  the  semi-circle  towards  the  wave  paddle.  The  elevation  of  the  turret 
was  adjustable  by  a rack  and  pinion  from  an  upper  (calibrate)  position 
above  the  water  surface  to  approximately  1 m beneath  the  surface;  and  it 
could  be  locked  in  any  position.  Elevation  control  was  provided  by  a ver- 
nier scale  on  the  rack. 

To  calibrate  the  probes,  the  strut  could  be  rotated  90°  and  the 
turret  locked  with  the  probe  tips  just  beneath  the  bottom  orifice  of  the 
calibrating  cylinder,  from  which  a water  jet  flowed  over  the  tips  at  a 
metered  rate.  Because  these  experiments  were  unique  and  very  high  meas- 
urement accuracy  was  desired,  calibration  procedures  and  error  limits  are 
described  in  detail  in  Appendix  A. 


3.4  Elevation  Measurements 


For  the  purpose  of  obtaining  wave  profiles  and  their  speed  of  travel 
(phase  velocity,  C) , measurements  of  surface  elevation  were  obtained  by  two 
complimentary  methods;  recording  wave  staffs  and  side-window  photography. 

Wave  staffs  consisted  of  120  cm  tapered  strips  of  circuit  board  stock 
exposed  edgewise  to  the  waves  and  bearing  a vertical  array  of  96  gold- 
plated  shorting  contacts  at  1 cm  intervals.  A logic  circuit  at  the  top  of 
each  staff  produced  an  analog  voltage  which  changed  by  125  mv  per  immersed 
contact  with  a response  rate  of  20  khz.  Since  meniscus  effects  were  limited 
to  a milimeter  or  so,  the  staff  resolution  was  1 cm,  except  in  breaking 
waves,  in  which  case  the  output  provided  only  an  upper  bound  to  surface 
elevation.  The  staffs  were  calibrated  during  data  taking  by  generating  a 
continuous  train  of  waves  of  maximum  height  and  examining  their  respective 
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contact  shorting  cycles  on  an  oscilloscope,  whose  dispaly  indicated  faulty 
contacts  or  incorrect  bridge  excitation  voltages  as  "missing  teeth". 

For  most  . c this  study,  three  staffs  were  mounted  on  the  rolling 
cart:  two  crosswise  of  the  channel  with  their  contact  edges  in  the  plane 

of  the  probe  tips  at  distances  of  9 cm  from  the  channel  sidewall  and  con- 
vergent barrier,  respectively;  and  one  staff  in  tandem  at  9 cm  from  the 
sidewall  and  displaced  variably  toward  the  wave  paddle  at  distances 
between  167  cm  and  280  cm.  The  two  transverse  staffs  provided  a means  of 
monitoring  cross-channel  variations  in  elevation  (almost  always  less  than 
2 cm).  The  tandem  staffs  were  generally  displaced  one-half  wavelength 
apart  at  the  operating  wave  frequency,  from  which  phase  velocity  could  be 
deduced  from  consecutive  arrival  times  of  a long  train  of  identical  waves. 
During  other  series  concerned  with  measurement  of  fluid  velocities,  as  many 
as  fou^  additional  staffs  were  employed,  both  on  the  cart  and  at  fixed 
locations  between  the  cart  and  the  wave  paddle. 

Although  limited  to  the  single  channel  viewing  window,  motion  picture 
and  single-frame  photography  provided  a means  of  testing  the  accuracy  of 
the  wavestaff  elevations,  as  well  as  demonstrating  that  a time  record  of 
surface  elevation  at  a fixed  point  is  practically  equivalent  to  a wave 
profile  in  space  centered  on  that  point  (see  5.2). 


3.5  Data  Acquisition 


All  analog  sensors  were  serially  multiplexed  at  the  rate  of  100  hz  per 
channel  through  an  A/D  converter  and  stored  on  magnetic  tape  in  the  IBM  1130 
computer.  At  a wave  phase  velocity  of  280  cm/sec,  this  data  rate  produced 
one  elevation  and  one  velocity  sample  for  each  2.8  cm  of  wave  travel  dis- 
tance. All  analog  signals,  including  that  from  the  paddle  arm  transducer 
were  also  routinely  monitored  on  a Statos  Mark  3 strip  chart  recorder, 
principally  as  a guide  to  sensor  function. 

In  order  to  minimize  redundant  data  storage,  data  acquisition  was 
limited  to  the  first  6-8  equilibrium  waves  of  each  series.  This  was 
accomplished  by  constructing  a digital  timer  having  an  adjustable  time 
delay  which  could  be  preset  in  0.01  second  intervals  from  zero  to  120  sec. 
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The  timer  was  started  by  a microswitc  >.  on  the  paddle  arm,  and  provided  the 
computer  with  a data  start  pulse  after  the  preset  time  delay.  Such  preci- 
sion time  control  provided  easy  identification  of  individual  waves  in  separ- 
ate data  series  made  at  different  velocity  probe  elevations,  and  permitted 
time-registering  of  data  samples  to  within  ±0.02  sec  during  the  construction 
of  velocity  field  contour  plots. 


4.  METHODOLOGY  AND  PROCEDURES 

These  experiments  are  roughly  divisible  into  two  categories  accord- 
ing to  procedures  and  objectives:  (a)  velocity  measurements  to  define  the 

flow  fields;  and  (b)  elevation  measurements  to  define  wave  profiles,  lead- 
ing to  estimates  of  potential  energy  flux  and  energy  loss  due  to  breaking. 
In  both  cases,  measurements  were  made  at  several  stations  before  breaking- 
-and,  in  the  latter,  at  several  post-breaking  stations--in  order  to  define 
the  rates  of  change  of  these  properties.  Similar  sets  of  measurements 
were  made  for  series  of  heavily  breaking  and  lightly  breaking  waves,  as 
determined  only  by  initial  wave  steepness.  Lastly,  in  some  cases  compara- 
tive measurements  were  made  at  three  different  wave  frequencies  (0.66, 

0.73  and  0.80  hz)  to  determine  whether  or  not  the  behavior  of  steep  waves 
obeys  Froude  similitude,  so  that  important  wave  parameters  are  scalable  in 
terms  of  frequency  only. 

Before  discussing  procedures,  the  concept  of  growth  rate  and  the 
distinction  between  light  and  heavy  breaking  need  clarification.  Prior  to 
the  measurements  described  herein,  a series  of  breaking  experiments  was 
performed  at  constant  wave  frequency.  Starting  with  a paddle  stroke  so 
small  that  no  wave  broke  in  the  test  section,  the  stroke  was  then  incre- 
mentally increased.  It  was  observed  that  breaking  first  occurred  farthest 
from  the  paddle  and  was  of  the  intense  (plunging)  type  (Fig.  5).  Progres- 
sively higher  strokes  resulted  in  moving  the  breaking  point  closer  to  the 
paddle  and  decreasing  the  breaking  intensity  (Fig.  6) . The  wave  height  at 
breaking  was  --if  anything--sl ightly  lower  for  barely  breaking  waves  than 
for  plunging  waves.  While  the  relative  constancy  of  breaking  height  can 
reasonably  be  attributed  to  a stability  limit  imposed  by  steepness  alone, 
the  variation  of  breaking  intensity  is  more  logically  associated  with 


Fig.  6a  Light  breaking  0.66  hz  wave  about  6 m Fig.  6b  Rear  view  of  heaving  breaking  waves, 

from  paddle. 


corresponding  variation  of  the  scaled  channel  convergence  rate,-L(3w/9x)/w, 
which  increases  inversely  with  distance  from  the  paddle.  If,  as  in  shallow 
water  (Van  Dorn,  1975),  energy  flux  is  conserved  prior  to  breaking,  then 
it  follows  that  breaking  intensity  .ts  directly  related  to  a wave's  growth 
history  within  the  convergence.  A fundamental  objective  of  the  experi- 
ments was  to  explore  the  nature  of  these  relationships. 


4.1  Velocity  Measurement  Procedures 


Having,  as  above,  experimentally  determined  the  paddle  strokes  (tape 
demodulator  voltages)  corresponding  to  light  and  heavy  breaking  at  the 
three  operating  frequencies,1  measurement  stations  were  established  where 
the  first  four  uniform  waves  in  each  series  had  achieved  the  following 
growth  stages: 

o Young  waves  of  relatively  low  steepness  (H/L  s 0.10),  at  the  station 
nearest  the  paddle  where  the  waveform  first  became  stable,  as  evi- 
denced by  the  absence  of  anomalous  peaks  in  the  Fourier  spectrum 
(App.  B). 

o Pre-b’-fc^king  waves  of  steepness  0.10  < H/L  < 0.12,  characterised  by 
markfc-i  profile  asymmetry,  and  where  videotape  photography  first 
exhibited  specific  evidence  of  crest  instability.  These  locations 
varied  by  a few  centimeters  from  wave  to  wave,  but  their  respective 
flow  fields  were  found  to  be  essentially  similar,  except  within  the 
upper  2-3  cm  beneath  their  crests. 

o Breaking  waves,  experimentally  determined  to  have  maximum  horizontal 
flow  velocities,  and  generally  characterized  by  maxima  cf  wave 
height  and  steepnesses  in  the  range  0.12  < H/L  < 0.14. 

These  definitions  are  rather  subjective,  particularly  as  regards  the 
distinction  between  light  and  heavy  breaking.  In  fact,  'breaking'  itself 
is  better  described  as  a continuous  development,  involving  increasing 
steepness,  height,  profile  asymmetry,  and  internal  velocities,  the  rate  and 
degree  of  which  arc  all  functions  of  frequency,  initial  steepness,  and 
position  along  the  convergence  (growth  rate).  Probably  the  most  useful  and 
easily  observable  definition  of  breaking  is  that  stage  of  development 
corresponding  to  maximum  wave  height. 


For  dimensional  consistency,  water  depths  were  adjusted  to  one-half  the 
smai  1-ampl itude  wavelength  (D  = g/4«f*) , despite  the  fact  that  Michc  (1944) 
has  shown  that  tinder  these  conditions  water  depth  has  an  insignificant 
influence  on  the  properties  of  steep  waves.  This  scaling,  however,  had  the 
advantage  that  the  paddle  stroke  amplitude  was  nearly  the  same  at  different 
frequencies  for  the  same  breaking  intensity. 
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Because  of  physical  limitations  of  the  channel,  it  was  not  possible 
to  produce  both  light  and  heavy  breaking  at  all  three  frequencies  at  all 
of  the  above  growth  stages.  For  this  reason,  as  well  as  because  of  the 
time  required  to  make  detailed  velocity  measurements,  efforts  were  con- 
centrated at  a frequency  of  f = 0.66  hz,  for  which  the  widest  range  of 
breaking  intensity  could  be  accommodated,  supplemented  by  additional  obser- 
vations of  heavily  breaking  0.73  and  0.80  hz  waves. 

Measurement  procedures  at  each  station  were  virtually  identical. 
Following  calibration,  the  probe  turret  was  initially  positioned  at  its 
lowest  possible  elevation  (88  cm  above  the  channel  bottom),  the  timer  was 
set  to  initiate  data  taking  with  the  arrival  of  the  first  uniform  wave,  and 
three  or  four  wave  series  were  recorded  to  establish  repeatability.  The 
turret  was  then  raised  by  4-5  cm,  and  the  data  sequence  repeated.  Elevation 
increments  were  reduced  as  the  wave  crest  elevation  was  approached  in  order 
to  better  resolve  high  velocity  gradients.  Because  the  probe-pair  accept- 
ance angle  was  limited  to  90°,  repeat  runs  were  conducted  at  the  lower 
turret  elevations  with  the  turret  rotated  by  ±45°,  so  as  to  resolve  velocity 
vectors  making  angles  greater  than  45°  to  the  horizontal.  During  all 
velocity  measurements,  surface  elevation  was  simultaneously  recorded  as 
previously  described. 


4.2  Elevation  Measurement  Procedures 

For  the  elevation  measurements,  ten  stations  along  the  21.4  m test 
section  were  selected  for  each  wave  frequency.  Station  spacing  was  densest 
near  the  breaking  points  of  pretest  waves.  These  stations  were  then 
occupied  in  pairs  by  deploying  two  to  three  cross-channel  wave  staffs  and 
one  or  two  tandem  staffs  at  a time,  and  moving  them  sequentially  between 
stations.  For  each  staff  configuration,  three  or  four  duplicate  wave 
series  were  generated,  starting  with  a paddle  amplitude  just  low  enough  to 
preclude  breaking,  and  increasing  the  stroke  incrementally  so  as  to  produce 
heavy,  intermediate,  and  light  breaking,  respectively.  The  frequency  was 
then  changed,  and  the  entire  schedule  repeated.  As  with  the  velocity 
measurements,  not  all  growth  stages  could  be  produced  for  all  frequencies, 
and  the  most  extensive  results  were  associated  with  0.66  hz  waves. 
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Because  of  station  overlap  between  the  velocity  and  elevation  measure- 
ments, a complete  list  of  station  identifiers,  paddle  distances  and  assoc- 
iated wave  heights  is  given  in  Appendix  D,  which  is  prefaced  with  remarks 
relevant  to  order  and  content. 


5.  RESULTS 

5.1  State  Parameters  for  Waves  of  Finite  Height 


The  work  of  Stokes  and  others  (Section  2.)  indicates  that  the  only 
independent  parameter  distinguishing  the  immediate  state  of  steady  waves  of 
finite  height  from  those  of  negligible  amplitude  is  wave  steepness,  H/L, 
whose  limiting  value,  H*/L*  - 1/7,  is  ordinarily  assumed  to  characterize 
breaking  inception.*  Since  the  present  experiments  dealt  with  unsteady 
waves,  whose  limiting  steepnesses  and  breaking  intensities  varied  signifi- 
cantly and  clearly  depended  upon  initial  steepness  and  growth  rate,  it  was 
early  apparent  that  some  additional  growth  parameter  was  necessary  to  define 
a wave's  immediate  rate  of  change. 

In  Nature,  wave  growth  is  dependent  upon  the  relative  amplitudes  and 
speeds  of  the  principal  Fourier  components  of  the  wave  field,  as  manifested 
by  the  'group'  phenomenon  for  a narrow  spectrum  (Section  6.2).  In  these 
experiments,  growth  was  produced  by  lateral  squeezing  of  individual  waves 
within  the  converging  channel.  The  scaled  channel  convergence  rate^  (in 
space)  is  given  by: 

G = -L*(8w/9x)/w(x) , (1) 

where  w(x)  is  the  channel  width  at  paddle  distance,  x.  Because  9w/8x  * 
-1/10  = constant,  and  w varied  inversely  with  x,  G was  a monotonically 
increasing  function  of  x as  a wave  progressed  down  the  channel,  and  fell 
within  the  range,  0.165  < G < 0.75,  for  all  three  experimental  frequencies. 

If  the  energy  flux  is  conserved  among  growing  waves,  then  wave  ener- 
getics should  be  determined  by  growth  rate  alone.  For  shoaling  periodic 

*L*  = I.2g/27tf2  js  thc  st0}<0S  limiting  wavelength  at  frequency,  f (Michel  1, 
1944). 

2 

Also  called  'growth  rate'  hereinafter. 
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waves.  Van  Dorn  (1975)  has  shown  that  potential  energy  flux  is  approximately 
conserved  right  up  to  breaking,  and  it  was  similarly  found  here  to  be  con- 
served in  deep  water.  A principal  result  of  these  experiments  is  to  show 
that  the  frequency,  steepness,  and  growth  rate  suffice  to  determine  wave 
kinematics  and  energetics,  and  to  examine  how  other  wave  properties  depend 
upon  these  parameters. 

Because  of  their  relevance  to  all  subsequent  discussion,  specific 
values  of  these  parameters  for  all  stations  considered  in  this  report  are 
listed  below  in  Tables  I and  II.  Both  tables  are  keyed  to-and  ordered  in 
the  same  manner  as  Table  D;  that  is,  for  each  wave  frequency,  wave  series 
A,  B,  C....are  ordered  in  terms  of  increasing  paddle  stroke  amplitude, 
measurement  stations  are  ordered  according  to  increasing  paddle  distance 
(and  growth  rate) . Other  tabulated  parameters  comprise  further  condensa- 
tions of  the  repetitive  series  averages  in  Table  D.  The  methods  of  averag- 
ing and  the  determination  of  error  estimates  are  given  in  App.  F.  The  95% 
confidence  limits  for  these  parameters  depend  slightly  upon  the  state  of 
wave  development,  and  are  expressed  below  as  percentages  of  the  tabulated 
values: 


Steepness 

Range 

<0.10 

0.10-0.14 

Post-breaking 


Steepness 

H/L 

±3% 

±4% 

+6% 


Potential 
Energy  Density 

±3% 

±3% 

±5% 


Potential 
Energy  Flux 

±5% 

±5% 

±6% 


Unless  otherwise  quajufied,  these  confidence  limits  are  reflected  in  future 
figures  by  error  bars  through  individual  data  points. 

In  both  tables,  frequency,  f,  is  that  of  paddle  motion,  later  shown 
to  be  that  of  the  principal  Fourier  component  (App.  B) . Steepness  was 
calculated  from  the  wave  height  averages  in  Table  D,  divided  by  the  Stokes 
Sth-order  wavelength  (Section  5.1.1).  Growth  rate  was  calculated  from  the 
local  channel  width,  w(x) , the  constant  convergence  rate,  8w/Bx  =-0.097 
(elsewhere  given,  nominally  as  1/10),  and  the  Stokes  limiting  wavelength, 

L*  = 1.2g/2irf  . Values  of  potential  energy  density  are  averages  from 
Table  D,  and  potential  energy  flux  was  obtained  by  multiplying  them  by 
the  Stokes  5th-order  group  velocity  (Section  5.6.2), 


i 
i 


tacJiSsasss^F??^^^  ; 


In  Table  II,  small  numbers  1 and  2 appearing  to  the  left  of  some 
station  numbers  refer  to  early  and  late  stages  of  ’breaking1  waves, 
respectively,  as  distinguished  by  initial  instability  and  terminal  plung- 
ing. For  heavy  breaking,  these  distinctions  tend  to  coincide. 


Table  I 

State  Parameters  for  Velocity  Measurement  Stations 


Wave  frequency  f = 0.66  hz,  Depth  D = 187  cm 


Breaker 

Station 

Paddle 

Channel 

Growth 

Wave 

Wave 

Type 

Distance 

Width 

Rate 

Height 

Steepne: 

x(cm) 

w(cm) 

G 

H(cm) 

H/L 

Light 

E4 

9S9 

155 

0.266 

52 

0.126 

Heavy 

B3 

857 

166 

0.250 

40 

0.100 

Heavy 

B9 

1771 

78 

0.535 

56 

0.132 

Heavy 

BIO 

1848 

71 

0.592 

58 

0.136 

Wave  frequency  f = 

0.73  hz,  Depth  D = 154 

cm 

Heavy 

FI 

163 

163 

0.208 

33 

0.103 

Heavy 

F2 

1801 

75 

0.454 

48 

0.137 

% TABLE  II 

s' 


STEEPNESS,  POTENTIAL  ENERGY,  POTENTIAL  ENERGY  FLUX 


I 

; 

t 


Lr’ 


l 


i 


Wave  frequency  f : 

Station  Paddle 
Distance 
x(cm) 

= 0.66  hz, 

Channel 

Width 

w(cm) 

Depth  D = 

Growth 

Rate 

G 

187  cm. 

Steep- 

ness 

TT/Li 

Potential 

Energy 

Density 

E^cm1) 

Potential 
Energy  Flux 

F(cn,J/sec ; 

A 1 

357. 

214.5 

0.1939 

0.0624 

33.5 

0.44341E  04 

2 

593. 

191.7 

0.2169 

0.0667 

41.1 

0.55390E  04 

3 

833* 

168.6 

0.2468 

0.0679 

40.5 

0.54706E  04 

A 

1069. 

145.8 

0.2855 

0.0724 

47.0 

0.64699E  04 

5 

1186. 

134.5 

0.3096 

0.0763 

50.5 

0.70545E  04 

6 

1379. 

115.9 

0.3595 

0.0839 

62.2 

O.09879E  04 

7 

1422. 

111.7 

0.3730 

0.0820 

57.9 

0.82977E  04 

8 

1615. 

93.1 

0.4479 

0.0906 

74*6 

0.11122E  05 

9 

1766. 

78.5 

0.5315 

0.0945 

80.2 

0.12187E  05 

10 

2002. 

55.8 

0.7504 

0.1111 

113.6 

0.19685E  05 

B 1 

430. 

207.5 

0*2005 

0.0886 

73.3 

0.10828E 

05 

2 

586. 

192.4 

0*2162 

0.0934 

82.0 

0.12386E 

05 

3 

857. 

166.2 

0.2503 

0.1000 

92.5 

0.14465E 

05 

4 

1118. 

141.1 

0.2951 

0.1056 

99.2 

0. 1597  IE 

05 

5 

1388. 

115.0 

0.3623 

0.1131 

117.3 

0.19691E 

05 

6 

1500. 

104.2 

0.4000 

0.1208 

132.0 

0.23190E 

05 

7 

1508. 

’Co. 4 

0.4030 

0.1182 

122.8 

0.21239E 

05 

8 

1577. 

96.8 

0.4309 

0.1227 

140.0 

0.24354E 

05 

9 

1771. 

78.0 

0.5348 

0.1317 

149.6 

0.28094E 

05 

a i 10 

1848. 

70.4 

0.5915 

0.1362 

170.7 

0.32967E 

05 

C 1 

357, 

21  ».5 

0.1939 

0.1126 

111.5 

0.18660E 

0. 

2 

593. 

1V1.7 

0.2169 

0.1149 

128.5 

0.21794E 

Ov 

3 

833. 

168.6 

0.2468 

0.1165 

115.5 

0.19769E 

05 

4 

1069. 

145.8 

0.2855 

0.1251 

150.8 

0.27179E 

05 

5 

1186. 

134.5 

0.3096 

0.1250 

143.7 

0.25880E 

05 

a 6 

1379. 

115.9 

0.3595 

0.1330 

162.7 

0.30P06E 

05 

» 7 

1422. 

111.7 

0.373G 

0.1302 

161.0 

0.29954E 

05 

8 

1615. 

93.1 

0.4479 

0.1310 

177.4 

0.33172E 

05 

9 

1766. 

78.5 

0.5315 

0.1299 

178.2 

0.33092E 

05 

10 

2002. 

55.8 

0.7504 

0.1173 

156.0 

0.26937E 

05 

0 1 

357. 

214.5 

0.1939 

0.1184 

128.5 

0.22243E 

05 

2 

593. 

191.7 

0.2169 

0.1212 

143.8 

0.25318E 

05 

3 

833. 

168.6 

0.2468 

0.1227 

128.5 

0.22S12E 

05 

a » 4 

1069. 

145.8 

0.2855 

0.1253 

148,9 

0.26957E 

05 

5 

1186. 

134.5 

0.3096 

0.1245 

143.5 

0.25763E 

05 

6 

1379. 

115.9 

0.3595 

7 

1422. 

111.7 

0.3730 

0.1239 

154.1 

0.27568E 

05 

8 

1615. 

93.1 

0.4479 

9 

1766. 

78.5 

0.5315 

0.1111 

122.0 

0.20246E 

05 

10 

2002. 

55.8 

0.7504 

0.1100 

118.0 

0.19463E 

05 

TABLE  II  (continued) 

Wave  frequency  f = 0.6o  Jr/.,  Depth  D = 187  cm. 

Station  Paddle  Channel  Growth  Steep- 

Distance  Width  Rato  ness 

x(cm)  w(cm)  C ll/L 


Potcntia  1 
Gnergy 
Density 
E (cm1) 


Potential 
Eneri'v  Pin-. 

F(cm3/sec) 


•"  ' W 

I - ? 


1 

586. 

192.4 

0.2162 

0.1213 

140.7 

0.24778E 

05 

2 

688. 

182.6 

0.2279 

0.1150 

105.0 

0.17821E 

05 

3 

857. 

166.2 

0.2503 

0.1235 

138.2 

0.24669E 

05 

4 

959. 

156.4 

0.2661 

0.1258 

146*1 

0.26434E 

05 

5 

1118. 

141.1 

0.2951 

0.1272 

146.5 

0.26746E 

05 

6 

1388. 

116.0 

0.3623 

0.1299 

152.9 

0*283846 

05 

1 

357. 

214.5 

0.1939 

0.1208 

132.7 

0.23304E 

05 

2 

593. 

191.7 

0.2169 

0.1233 

145.7 

0.25977E 

05 

3 

833. 

168.6 

0.2468 

0.1227 

132.2 

0.23478E 

05 

4 

1069. 

145.8 

0.2855 

0.1243 

145.6 

0.26123E 

05 

5 

1186. 

134.5 

0.3096 

0.1198 

139.5 

0.24352E 

05 

6 

1379. 

115.9 

0.3595 

0.1222 

133.0 

0.23545E 

05 

7 

1422. 

111.7 

0.3730 

0.1118 

123.7 

0.20614E 

05 

8 

1615. 

93.1 

0.4479 

0.1154 

124.0 

0.21096E 

05 

9 

1766. 

78.5 

0.5315 

0.1106 

114.2 

0.18907E 

05 

10 

2002. 

55.8 

0.7504 

0.1038 

110.9 

0.18177E 

05 

1 

35  7. 

214.5 

0.1939 

0.1240 

140.0 

0.25064E 

05 

2 

593. 

191.7 

0.2169 

0.1241 

149.4 

0.26764E 

05 

3 

833. 

168.6 

0.2468 

0.1096 

115.2 

0.18966E 

05 

4 

1069. 

145.8 

0.2855 

0.1124 

122.3 

0.20456E 

06 

5 

1186. 

134.5 

0.3096 

0.1231 

141.7 

0.25235E 

C: 

6 

1379. 

115.9 

0.3595 

0.1004 

96.5 

0.15104E 

05 

7 

1422. 

111.7 

0.3730 

0,1085 

1 1 5 , <* 

0.18874E 

03 

8 

1615. 

93.1 

0.4479 

0.1040 

100.6 

0.16064E 

05 

9 

1766. 

78.5 

0.5315 

0.1071 

106.2 

0.17242E 

05 

10 

2002. 

55.8 

0.7504 

0.1148 

125.4 

0.21269E 

05 

TABLE  II  (continued) 


Wave  frequency  f = 0.73  hz,  Depth  D = 154  cm. 


Station 

Paddle 

Channel 

Growth 

Steep- 

Potential 

Potential 

Distance 

Width 

Rate 

ness 

Energy 

Energy  Flux 

x(cm) 

w(cm) 

G 

H/L 

Density 

Ep(cmJ) 

F(cm3/sec) 

A 

1 

517. 

199.1 

0.1708 

0.0601 

20.7 

0.24634E 

04 

2 

753. 

176.3 

0*1929 

0.0641 

23.8 

0.28756E 

04 

3 

914. 

160*7 

0.2116 

0.0714 

29.0 

0.35890E 

04 

4 

1094. 

143.4 

0.2373 

0.0736 

31.0 

0.38706E 

04 

5 

1180. 

135.1 

0.2520 

0.0772 

34.7 

0.44046E 

04 

6 

1330. 

12C.6 

0.2823 

0.0800 

36.7 

0.47218E 

04 

7 

1428. 

111.1 

0.3064 

0.0800 

38.7 

0.49725E 

04 

8 

1561. 

98.3 

0.3466 

9 

1664. 

88.4 

0.3858 

0.0853 

44.3 

0.58289E 

04 

10 

1797. 

75.5 

0.4518 

0.0888 

47.9 

0.64063E 

04 

B 

1 

517. 

199.1 

0.1708 

0.1029 

63.7 

0.91453E 

04 

2 

753. 

176.3 

0.1929 

0.1085 

79.9 

0. 1 1824E 

05 

3 

914. 

160.7 

0*2116 

0.1152 

75.7 

0.11637E 

05 

4 

1094. 

143.4 

0.2373 

0.1234 

98.0 

0.15797E 

05 

5 

1180. 

135.1 

0.2520 

0.1273 

99.5 

0.16428E 

05 

6 

1330. 

120.6 

0.2823 

0.1268 

99.5 

0.16391E 

05 

7 

1428. 

111.1 

0.3064 

0.1217 

93.5 

0.14917E 

05 

8 

1561. 

98.3 

0.3466 

2 

9 

1664. 

88.4 

0.3858 

0.1278 

108.1 

0.17916E 

05 

10 

1797. 

75.5 

0.4518 

0.1239 

107,2 

0.17351E 

05 

c 

1 

517. 

199.1 

0.1708 

0.1042 

69.2 

0.10002E 

05 

2 

753. 

176.3 

0.1929 

0.1131 

83.6 

0.12695E 

05 

3 

914. 

160.7 

0.2116 

0.1152 

78.2 

0.12021E 

05 

2 

4 

1094. 

143.4 

0.2373 

0.1300 

111.2 

0.18683E 

05 

5 

1180. 

135.1 

0.2520 

C. 1 295 

86.0 

0.14394E 

05 

6 

1330. 

120.6 

0.2823 

0.1284 

102.4 

0.17026E 

05 

7 

1428. 

111.1 

0.3C64 

0.1262 

98,5 

0.16152E 

05 

8 

1561. 

98.3 

0.3466 

0.1295 

109.5 

0.18327E 

05 

9 

1664. 

88.4 

0.3858 

0.1245 

100.6 

0.16345E 

05 

10 

1797. 

75.5 

0.4518 

0.1202 

99.1 

0.15672E 

05 

te- 
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TABLC  II  (continued) 

Wave  frequency  f = 0.73  hz,  Depth  D = 154  on. 

Station  Raddle  Channel 
Distance  Width 
x(cm)  w(cm) 


Growth 

Rate 

G 


Steep* 

nness 

H/L 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


517. 

753. 

914. 

1094. 

1180. 

133C. 

1428. 

1561. 

1664. 

1797. 


517. 

753. 

914. 

1094. 

1180. 

1330. 

1428. 

1561. 

1664. 

1797. 


Potential 

Energy 

Density 

En(cmJ) 


199.1 

0.1708 

0.1116 

79.0 

176.3 

0.1929 

0.1242 

93.9 

160.7 

0.2116 

0.1158 

84.2 

143.4 

0.2373 

0.1278 

107,5 

135.1 

0.2520 

0.1338 

116.7 

120.6 

0.2823 

0.1306 

112.2 

111.1 

0.3064 

0.1300 

111.0 

98.3 

0.3466 

0.1170 

95.0 

88.4 

0.3858 

0.1187 

93.1 

75.5 

0.4518 

o.im 

82.1 

199.1 

0.1708 

0.1188 

82.5 

176.3 

0.1929 

0.1231 

92.5 

160.7 

0.2116 

0.1158 

80.5 

143.4 

0.2373 

0.1234 

99.5 

135.1 

0.2520 

0.1256 

109.5 

120.6 

0.2823 

0.1186 

92.4 

111.1 

0.3064 

0.113* 

80.5 

98.3 

0.3466 

0.1122 

83.0 

88.4 

0.3858 

0,1101 

77,7 

75.5 

0.4518 

0,1093 

75.8 

-~'y  v' 335T *v  ^ 3*5^2?®!^! 


Potential 
Energy  Flux 

F(cm3/6ec) 


0.11888E  05 
0.15213E  05 
0.12987E  05 
0.17809E  05 
0.20078E  05 
0.18917E  05 
0.18641E  05 
0.14745E  05 
0.14607E  05 
0.12330E  05 


0.12938E  05 
0.14893E  05 
0.12409E  05 
0.16039E  05 
0.17894E  05 
0.14490E  05 
0.12239F  05 
0.12533E  05 
0.11596E  05 
0.U256E  05 


m 


Kb.-. 
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TABLE  II  (continued) 


Wave  frequency  f = 0.80  hz.  Depth  D = 132  cm. 


Station 

Paddle 

Distance 

x(cm) 

Channel 

Width 

w(cm) 

Growth 

Rate 

r. 

Steep- 

ness 

H/L 

Potential 

Energy 

Density 

Ep(cmJ) 

Potential 
Energy  Flux 

F(cm3  /sec) 

A 

1 

804* 

171.4 

0.1652 

0.0675 

18.7 

0.20864E 

04 

2 

1040* 

148.6 

0.1906 

0.0702 

20.4 

0.23020E 

04 

3 

1056. 

147.0 

0.1927 

0.0738 

22.7 

0.25*>41E 

04 

4 

1292. 

124.3 

0.2281 

0.0797 

29.3 

0.34325E 

04 

5 

1532. 

101.1 

0.2806 

0.0857 

32.0 

0.38440E 

04 

6 

1675. 

87.3 

0.3252 

0.0915 

35.5 

0.43825E 

04 

7 

1768. 

78.3 

0.3626 

0.0958 

39.9 

0.50338E 

04 

8 

1911. 

64.5 

0.4407 

0.1042 

46.8 

0.61802E 

04 

B 

1 

804. 

171.4 

0.1652 

0.1070 

49.5 

0.66252E 

04 

2 

1040. 

148.6 

0.1906 

0.1072 

50.2 

0.67276E 

04 

3 

1056. 

147.0 

0.1927 

0.1085 

52.5 

0.70847E 

04 

4 

1292. 

124.3 

0.2281 

0.1204 

66.0 

0.95468E 

04 

5 

1532. 

101.1 

0.2806 

0.1327 

79.7 

0.12428E 

05 

i 

6 

1675. 

87.3 

0.3252 

0.1307 

76.2 

0.11738E 

05 

a 

7 

1768. 

78.3 

0.3626 

0.1334 

77.7 

0.12175E 

05 

8 

1911. 

64.5 

0.4407 

0.1305 

79.3 

0.12188E 

05 

C 

1 

804. 

171.4 

0.1652 

0.1144 

59.0 

0.82289E 

04 

2 

1040. 

148.6 

0.1906 

0.1128 

56.2 

0.77817E 

04 

3 

1056. 

147.0 

0.1927 

0.1136 

55.2 

0.76741E 

04 

4 

1292. 

124.3 

0.2281 

0.1281 

77.5 

0.11746E 

05 

2 1 

5 

1532. 

101*1 

0.2806 

0.1365 

85.7 

0.13692E 

05 

6 

1675, 

87.3 

0.3252 

0.1307 

73.0 

0.11238E 

05 
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5.1.1  Wave  Propagation  Speed 

The  determination  of  wavelength,  steepness,  and  energy  flux,  as  well 
as  the  conversion  of  time-  to  space-scales,  all  require  a suitable  defini- 
tion of  the  instantaneous  speed  of  moving  waves.  For  steady,  irrotational 
waves  of  finite  height,  the  Stokes  surface  profiles  to  any  order  are 
symmetric  about  their  crests.  All  points  of  constant  phase  (crests  or 
troughs,  say)  move  at  the  same  speed  (phase  velocity,  C) , experimentally 
determinable  by  observing  the  transit  times  of  individual  phases  over  a 
known  distance.  Alternatively,  where  the  dispersion  relation,  C = f(H,T), 
is  known  analytically,  C can  be  determined  from  measurements  of  wave  height, 
H,  and  period,  T.  The  latter  is  ordinarily  taken  as  the  reciprocal  of  the 
paddle  frequency,  f.  For  Stokes  waves  of  3rd  and  5th  order,  respectively, 
phase  velocity  is  given  implicitly  by: 


• fe  » ♦ <tf >2> 


= ^ — [1 
2trf  11 


* CIS)2  * c-5S*)4l 

LS  L5 


(Levi-Civita,  1925)  (2a) 

(Beach  Erosion  Board,  1941)  (2b) 


In  the  present  experiments,  however,  all  waves  withir.  the  convergent 
channel  were  unsteady,  and  their  profiles  became  increasingly  asymmetric 
as  their  height  increased  towards  breaking.  The  degree  to  which  their 
instantaneous  speeds  fit  the  above  steady-state  models  was  investigated  by 
two  independent  methods: 

o crest  speeds  were  determined  by  timing  the  arrivals  of  corresponding 
elevation  maxima  at  two  tandem  wave  staffs  less  than  a wavelength 
apart; 

0 phase  speeds,  defined  as  the  velocity  of  wave  spectral  components 
associated  with  the  paddle  frequency,  f,  were  calculated  as  follows. 

Simultaneous,  25-second,  surface  elevation  series  from  two  tandem  wave 
staffs  were  Fourier-analyzed,  and  their  cross-spectra  (Jenkins  and  Watts, 
1968)  were  averaged  over  triplicate  wave  series.  Phase  speed  was  then 
determined  from  the  relation,  C = 2irfAx/$,  where  Ax  was  the  staff  spacing, 
and  $ the  phase  of  the  cross  spectrum  at  frequency,  f.  This  latter  method 
had  the  advantage  of  demonstrating  that  the  paddle-generated  waves  were 
essentially  free  from  energy  peaks  at  frequencies  other  than  f and  its 
higher  harmonics  (App.  B) . 
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Figure  7 compares  representative  values  of  phase  and  crest  speeds  so 
obtained*  for  thirty  different  wave  heights,  divided  among  three  paddle 
frequencies,  with  curves  computed  from  equations  2a  and  2b,  from  which  the 
following  conclusions  were  drawn 

o The  scatter  of  phase  speed  was  somewhat  greater  than  that  for  crest 
speeds,  although  nowhere  greater  than  ±3%  about  the  median  speed 
for  a particular  wave  height.  This  result  is  somewhat  surprising. 
One  light  expect  that  phase  speed  would  be  more  representative  of 
wave  propagation  because  spectra  indicate  that  over  70%  of  wave 
energy  is  concentrated  in  a band  0.04  hz  wide,  centered  on  the 
paddle  frequency,  whereas  the  crest  is  a point  on  a wave  whose 
position  can  fluctuate  due  to  small  disturbances. 

o The  locus  of  crest  speeds  is  closer  to  the  Sth-order  Stokes  curve 
than  to  the  3rd-order  curve,  although  these  curves  differ  by  only 
1%  at  the  Stokes  (Michell)  steepness  limit,  H*/L*  = 0.142.  Owing 
to  the  nature  <f  the  power  expansions  on  which  these  curves  are 
based,  neither  passes  precisely  through  the  limit  phase  speed  at 
limiting  steepness.  However,  the  locus  of  phase  speeds  (not  for- 
mally determined)  does  appear  to  pass  through  this  point,  and  lies 
consistently  below  that  for  crest  speeds.  This  is  consistent  with 
the  premise  that  crest  speeds  should  progressively  exceed  phase 
speeds  in  growing  waves,  and  that  the  latter  should  better  agree 
with  a proper  theory. 

However,  because  the  overall  error  is  small  in  any  case,  and  because 
the  5th  order  theory  gives  the  best  approximation  to  observed  crest  speeds, 
it  is  used  hereinafter  for  tne  calculation  of  phase  velocity  and  wavelength 
up  to  breaking.  It  is  also  used,  where  necessary,  after  breaking,  although 
its  applicability  is  not  established  because  of  the  uncertainty  of  deter- 
mining crest  arrivals. 

To  determine  whether  the  scatter  of  phase  speeds  in  Fig.  7 might  be 
due,  in  part,  to  finite  growth  rate,  the  influence  of  frequency  and  steep- 
ness was  removed  by  normalization  with  respect  to  the  5th-order  phase 
velocity,  C,.,  and  the  results  plotted  vs  growth  rate  in  Fig.  8.  The  lack 
of  any  apparent  trend  implies  a negligible  influence  of  G,  compared  to  that 
of  frequency  or  steepness. 


*The  abscissa  for  each  datum  is  an  average  over  the  first  four  waves  of 
triplicated  series.  Crest  speed  ordinates  are  averages  over  about  45 
consecutive  waves.  The  two  representative  error  bars  denote  84%  confidence 
limits  from  an  assumed  Gaussian  distribution  of  variances  in  arrival  times. 
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Fig.  8 Influence  of  growth  rate  upon  wave  speed. 
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5.1.2  Steepness-Growth  Rate  Diagram 

It  was  earlier  remarked  (Section  4)  that,  while  breaking  appeared  to 
be  governed  by  some  limiting  steepness,  breaking  intensity  was  qualitatively 
linked  to  growth  rate.  These  relationships  can  be  quantified  by  plotting 
H/L  vs  G for  individual  waves,  since  any  point  on  such  a diagrar:  defines 
the  immediate  state  of  a wave,  and  a sequence  of  such  points  determines  a 
trajectory  that  gives  the  wave's  steepness  history  as  it  grows  toward 
breaking. 


a 


r 


Figure  9 shows  two  representative  histories  for  light  and  heavily 
breaking  0.66  hz  waves,  corresponding  to  data  series  E and  B in  Table  II, 
respectively. * Error  bars  indicate  the  steepness  uncertainty  attendant  to 
data  averaging.  Each  trajectory  commences  at  the  left  at  some  initial 
steepness  determined  by  paddle  stroke  amplitude,  and  slopes  upward  to  the 
right  as  the  wave(s)  converge  towards  breaking,  as  identified  by  the 
numbers  1 or  2 previously  defined.  In  this  representation,  light  breaking 
is  evidently  distinguished  from  heavy  breaking  by  a higher,  initial  steep- 
ness, and  by  a lower  growth  rate  at  breaking.  Visually,  the  difference  in 
breaking  intensity  is  obvious  (Figs.  5 C 6),  and  is  later  shown  to  depend 
directly  upon  growth  rate  (Section  5.7.1).  The  fact  that  both  curves  appear 
not  to  have  reached  maximum  steepness  at  breaking  may,  in  part,  be  due  to 
the  difficulty  of  determining  steepness  accurately  after  breaking.  But  it 
also  can  be  qualitatively  argued  that  breaking  within  a uniform  convergence 
represents  some  sort  of  equilibrium  between  energy  supplied  by  lateral 
squeezing  and  that  lost  by  breaking,  and  that  wave  height-  and  hence 
steepness-  may  actually  increase  slightly  beyond  the  breaking  point  (see 
Section  5.  7. 2) . 


5.1.3  Breaking  Steepness 

The  two  curves  in  Fig.  9 arc  only  representative  examples  of  an 
infinite  set,  whose  initial  steepnesses  range  from  zero  to  the  theoretical 
limiting  value,  0.142,  and  whose  maxima  are  bounded  by  this  limit,  or  by 

^Letter-number  symbols  in  this  figure  identify  wave  states  later  disc ’.ssed 
in  other  context. 
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Typical  steepness-growth  rate  trajectories  for  light  and  heavy  breaking  0.66  hz  waves 
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some  empirical  lower  bound.  To  better  define  this  latter  bound.  Fig.  10 
shows  the  maximum  (breaking)  steepnesses  for  all  wave  series  in  Table  II 
plotted  against  growth  rate.  Despite  the  large  error  bars,  which  are  the 
result  of  expanding  the  steepness  scale  within  the  breaking  range,  there 
appears  to  be  a significant  tendency  for  breaking  steepness  to  increase 
with  growth  rate.  This  is  opposite  to  the  expectation  that  those  waves 
which  grow  slowest  should  achieve  maximum  steepness-possibly  approaching 
the  theoretical  limit  (0.142)  for  infinitely  slow  growth.  The  converse 
result  implies  some  steepness  moderating  factor  near  breaking  that  is  only 
transcended  by  rapid  growth.  A possible  explanation  is  that  steep  waves 
are  weakly  dissipative  before  breaking  because  of  strong  velocity  gradients 
beneath  their  crests  (Section  5.4.2),  and  would  never  achieve  limiting 
steady- state  steepness  however  small  their  growth  rates.  If  breaking 
energy  loss  is  a monotonical ly  increasing  function  of  steepness,  equili- 
brium breaking  steepness  should  increase  with  growth  rate,  as  observed. 

This  point  is  further  elaborated  in  Section  5.6. 

5. 2 Stationarity  in  Growing  Waves 

Most  of  the  operational  data  in  these  experiments  were  obtained  from 
time-series  records  made  at  fixed  stations.  A mathematically  more-tractiblc 
representation  of  wave  properties,  that  can  be  compared  to  steady  wave 
theory,  is  one  referred  to  as  stationary  spatial  coordinate  system.  A meas- 
ure of  stationarity  is  the  accuracy  with  which  a stationary  time  series  can 
be  equated  to  synoptic  observations  at  an  instant.  This  question  was 
examined  here  by  comparing  staff  records  of  surface  elevation  to  simultan- 
eous motion  pictures  of  steep  (H/I,  = 0.136),  breaking,  0.66  hz  waves  at 
high  growth  rate  (G  = 0.54).  These  conditions  represent  the  most  rapid 
changes  observed,  and  thus  the  most  critical  test  of  stationarity. 

Figure  11  compares  three  different  superimposed  profiles  of  the  same 
wave.  The  solid  line  is  the  true  space  profile,  reconstructed  from  a 
single- frame  photograph  at  the  instant  when  the  moving  wave  was  most  nearly 
centered  on  the  wave  staff.  The  abscissa  of  the  frame-by-frame  and  wave 
staff  records  were  obtained  by  mult, plying  corresponding  phase  arrival  time 
increments  by  the  Stokes  I -oj^r  phase  velocity,  Cj.  (Oq.  2b).  Fxcept  on 
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the  overhanging  crest  face,  where  the  wave  staff  indicated  a near- vertical 
rise,  all  three  profiles  agree  within  1.5  cm  (2%  of  CgAt) . This  comparison 
also  served  to  verify  the  accuracy  of  wave  staff  response  and  resolution. 


5.3  Scaling  and  Similitude 

In  order  to  compare  or  contrast  flow  fields  and  surface  profiles  of 
different  waves  it  is  necessary  to  know  if  they  are  dynamically  similar. 
This  requires  the  following  conditions  be  satisfied  (Duncan,  et  al . , 1960, 
pp.  178-185) 

o Corresponding  space-time  positions  are  related  by  constant  linear 
scale  factors,  Xj/d^  = *2^2’  yl^l  = ^2^2*  Zl^dl  = z2^2'  and 

Vl  = t2e2* 

o All  relevant  physical  properties  and  locations  of  the  boundaries 
must  be  distributed  at  corresponding  positions. 

o Forces  or  scalar  quantities  at  corresponding  positions  can  be 
related  by  constant  linear  scale  factors. 

For  all  frequencies  in  these  experiments  the  dominant  forces  were 

gravitational  and  inertial,  and  it  was  convenient  to  let  d = L*  and  e = f 

be  linear  space  and  time  scales  respectively.1  The  second  condition 

requires  w(x)/L*  = constant,  9w/3x  = constant,  and  D/L*  = constant,  as  was 

? 

the  case  in  this  study."  The  third  condicion  is  a corollary  of  the  con- 

2 

stancy  of  Froude  number,  C /L*g  (ibid,  p.  182).  It  was  shown  (Section 

5.1.1)  that  C,  the  wave  propagation  speed,  can  be  expressed  as  the  product 

of  a function  of  steepness,  F(H/L) , and  (g/2Tif) , where  L = C/f.  Recalling 
2 

that  L*  = 1.2g/27if  , it  is  clear  that  the  Froude  number  is  proportional  to 
F(H/L),  so  that  steepness,  H/L,  must  also  be  constant  for  dynamical  simil- 
itude. For  the  frequencies  used  in  this  study,  the  nominal  values  of  the 
constant  scale  factors  are: 


^his  choice  of  scaling  parameters  had  the  advantage  that  when  the  waves 
were  breaking,  C/C*  = L/L*  = 1. 

2 

This  is  equivalent  to  growth  rate,  G = -L*(3w/9x) /w(x)  = constant. 
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Frequency 

f 

(hz) 

Length 

L* 

(cm) 

Speed 

c* 

(cm/sec) 

Energy 

gpLfxl0-8 

(ergs/cnr) 

I 

i 

\ 

4 

0.66 

430 

284 

1.82 

0.73 

352 

257 

1.22 

/ 

0.80 

293 

234 

0.84 

In  summary,  with  space  and  time  scale  factors,  L*  and  f,  respectively 
it  is  sufficient  to  maintain  constant  growth  rate,-L*8w/8x/w(x) , and  steep- 
ness, H/L,  in  order  to  achieve  dynamical  similitude. 


5 . 4 Results  of  Flow  Velocity  Measurements 


The  measurement  of  internal  flow  velocities  in  growing  waves  was  one 
of  the  principal  objectives  of  these  experiments,  not  only  because  they 
comprise  the  most  sensitive  test  of  any  theory,  but  also  because  they  per- 
mit some  insight  of  the  breaking  process,  for  which  there  presently  is  no 
adequate  theory.  As  described  in  Section  3.3,  flow  measurements  were  con- 
fined to  velocities  having  components  in  the  direction  of  wave  propagation, 
to  accuracies  of  order: 

Range  Magnitude  Direction 

cm/sec  cm/sec 

30<U<65  +6  4° 

65<U<550  3%  2° 

Flow  velocities  are  discussed  here  in  two  different  contexts.  In  section 
5.4.1,  profiles  of  maximum  horizontal  velocity  vs  elevation  beneath  indivi- 
dual wave  crests  are  compared  to  those  computed  for  steady  waves  of  corres- 
ponding steepness  from  Stokes  higher-order  theory.  This  representation 
serves  to  maximize  the  discrepancy  between  theory  and  observation  as 
unsteady  waves  grow  towards  breaking,  and  to  emphasize  the  influence  of 
growth  rate  on  breaking  intensity.  Section  5.4.2  represents  the  internal 
velocity  fields  in  growing  waves  in  the  form  of  isolines  (contours)  of 
constant  velocity  in  space-time  coordinates,  upon  which  are  superimposed 
direction  arrows  and  measured  wave  profiles.  Such  plots  graphically 
illustrate  how  increasing  steepness  ar.d  profile  asymmetry  are  reflected 
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internally  by  distortion  of  the  velocity  field  in  a manner  not  predictable 
for  steady  waves. 


5 . > Profiles  of  Horizontal  Velocity 


Although  detailed  measurements  of  internal  flow  velocities  were  con- 
ducted at  over  forty  stations,  six  exemplary  proxies  of  maximum  horizontal 
velocity  suffice  to  illustrate  without  redundancy  the  principal  features  of 
wave  development  as  functions  of  frequency,  steepness,  and  growth  rate. 

These  six  stations  are  identified  by  letter-number  symbols  in  the  steepness- 
growth  rate  diagram  (Fig.  9),  and  relevant  wave  parameters  are  listed  in 
Table  I. 

If  steady  wave  theory  has  any  application  to  unsteady  waves,  it  should 
best  apply  to  (young)  waves  of  low  steepness  and  growth  rate.  For  these 
conditions.  Fig.  12  shows  plots  of  the  vertical  distribution  of  maximum 
horizontal  velocity  for  the  first  four  uniform  waves  in  heavy-breaking 
series  B-3  and  F-l,  at  frequencies,  0.66  and  0.73  hz,  respectively.  Also 
plotted  are  other  data  for  waves  of  similar  steepness,  but  higher  frequency 
(1.15  hz) , reported  by  Morison  and  Crooke  (1953). 1 The  solid  curves  in 
this  figure  were  computed  from  Stokes  5th-order  theory  (Skjelbreia  and 
Hendrickson,  1960),  using  appropriate  water  depths  and  steepnesses.  In 
this-  and  the  two  succeeding-  figures,  measurement  elevations  and  flow 
velocities  have  been  normalized  to  the  Stokes  limiting  wavelengths  and 
phase  velocities,  respectively.  For  constant  steepness,  this  Froude  scal- 
ing should  eliminate  the  influence  of  frequency,  except  in  so  far  as  it 
affects  growth  rate;  that  is,  profile  differences  should  be  functions  only 
of  growth  rate. 

As  an  engineering  approximation,  the  computed  curves  might  be  regarded 
as  a reasonable  upper  bound  to  the  observed  data,  but  there  are  consistent 
differences  that  exceed  the  above-stated  measurement  accuracy.  The  fit  is 
acceptable  near  the  crests  and  well  below  trough  level,  but  the  velocities 
observed  here  are  uniformly  low  in  the  mid-range  by  about  25-30%  (12-15 
cm/sec),  and  their  vertical  distribution  is  more  hyperbolic  than  exponential, 

^orison  and  Crooke  computed  velocities  from  photographs  of  neutral-density 
oil  droplets.  Because  of  slightly  lower  steepness  (0.086  vs  0.10),  their 
normalized  wave  height  range  was  about  14%  smaller,  as  indicated  by  the 
dashed  lines  in  Fig.  11. 
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Fig.  12  Maximum  horizontal  flow  velocities  within  low  steepness  waves. 
Open  circles  represent  data  from  Mori  son  and  Crooke  (1953). 
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as  predicted  by  the  Stokes  curves.  As  shown  below,  this  tendency  increases 
with  steepness.  Secondly,  the  growth  rate  is  about  25%  higher  for  the  0.66 
hz  waves  than  for  the  0.73  hz  waves,  whereas  the  measured  velocities  are 
slightly  lower-particularly  in  mid-range.  This  result  is  opposite  to  that 
observed  at  higher  steepnesses,  but  the  difference  is  within  measurement 
error,  and  could  have  resulted  from  the  somewhat  arbitrary  choice  of  the 
probe  response  exponent,  n , in  the  low-velocity  range  (App.  A).  Similar 
differences  between  the  present  data  and  those  of  Morison  and  Crooke  are 
inexplicable,  except  to  note  that  they  are  opposite  to  what  would  be 
expected  for  waves  of  smaller  steepness  and  (presumably)  zero  growth  rate. 

But  these  comparisons  are  pushing  the  threshold  accuracy  of  the 
experiments,  and  the  principal  conclusion  is  that  flow  velocity  in  waves 
of  steepness  <0.1  are  predicted  acceptably  by  Stokes  5th-order  theory. 

Figure  13  is  a similar  plot  for  pre-breaking  waves  of  the  same  fre- 
quency and  nearly  the  same  steepness  (series  E-4  and  B-9,  Fig.  9 and 
Table  I).  The  comparison  here  is  the  difference  between  the  velocity  pro- 
files within  lightly-  and  heavily -breaking  waves  near  terminal  growth,  as 
manifested  by  a large  difference  in  growth  rate,  G.  Stokes  5th-order 
curves  for  both  steepnesses  are  given,  together  with  that  for  limiting 
steepness  (H*/L*  = 0.142),  computed  from  expansion  coefficients  given  by 
Michel  1 (1895). 

Again,  both  data  distributions  are  in  substantial  agreement  below 
trough  level,  but  increasingly  diverge  upward.  Near-crest  velocities 
within  waves  destined  to  break  heavily  (high  growth  rate)  exceed  those  for 
light  breaking  by  as  much  as  40%.  Neither  distribution  much  resembles  any 
of  the  computed  curves,  all  velocities  being  generally  lower  until  just 
beneath  crest  height,  and  then  increasing  hyperbolical ly  so  as  to  exceed 
even  those  predicted  by  the  limiting  wave  profile.  Thus,  for  waves  of 
similar  near-breaking  steepness,  the  influence  of  increasing  growth  rate  is 
to  markedly  increase  the  vertical  gradient  of  maximum  horizontal  velocity 
and,  as  later  shown,  to  correspondingly  increase  breaking  intensity. 

Lastly,  Fig.  14  compares  observed  profiles  for  two  heavily  breaking 
waves  of  the  same  steepness,  but  different  frequency  (series  B-10  and  F-2, 
respectively,  Fig.  9),  to  the  Stokes  limiting  wave  profile.  The  close 
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Fig.  13  Comparison  of  maximum  horizontal  flow  vcloc  ties  at  low  (o)  and 
high  (x)  growth  rates  near  breaking. 
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data  correspondence  between  frequencies  again  supports  Froude  similitude— 
even  at  breaking.  Crest  flow  velocities  slightly  exceed  the  limiting  phase 
velocity,  C*,  although  theoretical  limiting  steepness  is  never  quite 
achieved.  If  the  25%  higher  growth  rate  of  0.66  hz  waves  has  any  influence 
on  velocity  distribution,  it  is  to  slightly  lower  the  elevation  at  which 
maximum  velocity  occurs.  This  feature  is  further  discussed  in  the  next 
section. 

Although  the  observed  velocity  distributions  deviate  markedly  from 
the  computed  limiting  Stokes  profile,  they  appear  to  correspond  quite 
closely  to  those  observed  by  Van  Dorn  (1975)  within  shallow-water  waves 
breaking  on  gentle  slopes.  Irrespective  of  frequency,  he  found  that  the 
ratio  of  maximum  horizontal  velocity,  U,  to  observed  phase  velocity,  C,  was 
given  by  the  hyperbolic  equation: 


U/C  = aZ/(b  - Z)  + UQ/C,  (3)  ^ 

where  Z « Y/Y  was  the  fractional  ordinate  (bottom  to  crest)  at  which  U was 

Iu 

measured,  a = 0.10  and  b = 1.125  are  arbitrary  constants,  and  UQ/C  * 0.20 
at  the  bottom.  In  Fig.  14,  U/C  * 0.20,  approximately,  at  trough  elevation, 
zt,  for  which  the  appropriate  normalization,  referred  to  still  water  level, 
requires  that  Z * (z-zt)/H.  The  heavy  curve  in  Fig.  13  was  computed 
from  Eq.  (3),  using  the  same  coefficients.  Although  its  extrapolation 
below  trough  level  is  obviously  unwarranted,  the  fit  above  this  level  sug- 
gests that  the  breaking  of  deep  and  shallow  waves  have  marked  dynamical 
similarities. 


I 


5.4.2  Internal  Velocity  Fields 

Although  over  600,000  independent  determinations  of  fluid  velocity 
and  direction  were  obtained  in  the  course  of  these  experiments,  a substan- 
tial fraction  comprised  redundant  repetitions  to  insure  accuracy.  It  was 
also  found  that,  while  discrete  data  were  highly  reproducible  at  all  points 
within  individual  waves  of  repetitive  and  quasi-uniform  series,  the  details 
of  the  velocity  distributions  near  the  crests  of  steep  waves  varied  just 
enough  from  wave  to  wave  sr  as  to  obscure  the  fine  structure  which  it  was 
desireable  to  resolve.  This  result  proscribed  near-crest,  wave-to-wave 
averaging  of  velocity  data  in  the  same  manner  as  employed  for  wave  height 
or  steepness.  The  fact  that  individual  flow  representations  of  the  type 
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described  below  looked  very  much  alike  from  wave  to  wave  argues  that  the 
above  differences  are  attributable  to  slight  phase  shifts,  rather  than  to 
timing  errors.  At  depths  greater  than  about  22  cm  beneath  steep  crests, 
flow  details  were  sufficiently  similar  so  as  to  permit  meaningful  averages. 
Accordingly,  the  few  examples  included  here  were  selected  from  a very 
large  body  of  data  as  typically  illustrative  of  the  instantaneous  velocity 
fields  at  principal  stages  of  wave  growth. 

All  velocity  fields  were  generated  from  computer  print-outs  in  which 
flow  velocity  magnitudes  were  displayed  as  numbers,  vertically  arrayed  in 
proportion  to  scaled  measurement  elevation,  and  horizontally  according  to 
synchronized  time  sequence.  Horizontal  number  spacing  was  scaled  according 
to  phase  velocity  so  that  the  array  corresponded  to  an  undistorted  x-z 
field  at  the  instant  of  crest  passage.  Flow  directions  were  indicated  by 
a computer- generated  arrow  between  each  pair  of  numbers.  Smoothed  contour 
lines  of  constant  velocity  were  then  drawn  by  eye  across  each  array.  For 
elevations  lower  than  22  cm  beneath  a wave  crest,  arrays  for  the  first  four 
uniform  waves  in  a series  were  superimposed  to  better  define  average  flow; 
above  this  level  only  that  array  having  the  highest  velocities  was  contour- 
ed. Each  plot  was  bounded  upwardly  by  the  corresponding  synchronized 
surface  elevation  record  from  staff  data.  Considering  that  each  array 
spanned  only  about  120°  of  wave  phase  and  comprised  upwards  of  3,000  data 
points,  it  is  felt  that  even  the  fairly-complicated  flow  detailing  portray- 
ed at  breaking  is  an  accurate  representation.  Accuracy  was  considered 
lowest  where  flow  directions  were  inclined  more  than  45°,  so  as  to  produce 
backflow  on  either  probe  tip,  although  the  directional  calibrations  were 
still  consistent  and  reproducible.  These  areas  are  indicated  by  shading  in 
the  following  figures. 

Figures  15,  16,  and  17  show  velocity  fields  so  constructed  for  young, 
pre-breaking,  and  breaking  stages  of  a heavily-breaking  0.66  hz  wave, 
corresponding  to  data  series  B-3,  B-9,  and  B-10,  respectively,  in  Fig.  9 
and  Table  I.  The  profile  of  the  young  wave  is  fairly  symmetric  about  its 
crest,  as  are  its  internal  flow  directions  and  velocity  contours.  The 
latter  are  widely  spaced  and  convex  upward.  As  shown  in  Fig.  12,  even 
maximum  horizontal  velocities  are  close  to  those  computed  from  Stokes  5th- 
order  theory,  and  a velocity  field,  so  computed,  would  undoubtedly  look 
very  similar  to  Fig.  15. 
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Figure  16  shows  the  same  wave  about  0.25  sec  (0.18  wavelengths) 
before  breaking.  Steepness  has  increased  by  32%,  and  the  surface  profile 
has  become  markedly  asymmetric.  The  higher  velocity  contours  are  now 
concave  upward  and  concentrated  toward  the  forward  crest  face.  All  con- 
tours are  tilted  in  the  direction  of  wave  advance,  and  only  the  flow 
directions  resemble  those  in  the  previous  illustration.  The  lowest  veloc- 


ity contour  shown  (80  cm/sec)  is  also  the  highest  shown  in  Fig.  15. 

Figure  17  shows  the  wave  at  the  point  of  plunging.  The  upper  veloc- 
ity contours  have  now  wrapped  around  a high-velocity  jet  issuing  from  the 
forward  face  just  beneath  the  crest  at  about  1.1  times  phase  velocity.^ 

The  crest  face  itself  is  approaching  verticality,  and  the  flow  directions 
are  strongly  skewed  toward  the  issuing  jet.  This  picture  is  in  great 
contrast  to  a similar  representation  for  a Stokes  limiting  wave  of  the 
same  height  and  frequency  (Fig.  18),  whose  flow  directions  and  velocity 
contours  were  computed  from  Michel l 's  expansion  coefficients.  More  than 
any  other  evidence,  this  comparison  demonstrates  the  inapplicability  of 
steady  wave  theo’  • to  unsteady  waves.  While  the  stream- function  method 
of  Dean  (1965)  was  developed  to  calculate  the  interior  potential  fields 
within  steady,  irrotational  waves  of  arbitrary  surface  profile,  Dean's 
method  must  fail  at  breaking,  where  the  free  surface  is  no  longer  a stream-  ' 
line,  and  the  interior  velocities  are  double-valued  along  the  vertical. 

While  these  details  are  of  obvious  significance  to  an  understanding  of  the 
breaking  process,  they  comprise  such  a small  fraction  of  the  total  field 
of  motion  as  to  be  of  only  secondary  importance  to  the  calculation  of 
forces  and  their  moments  for  engineering  purposes.  It  remains  to  be  shown 
whether  Dean's  method  otherwise  predicts  the  essential  velocity  structure. 

The  distinction  between  the  interior  flow  fields  of  heavy  vs  light- 
breaking waves  can  be  seen  by  comparing  Figs.  16  and  19.  Both  waves  have 
about  the  same  height,  steepness,  surface  profile,  and  are  represented 
about  0.18  wavelength  before  breaking.  However  both  growth  rate  and  flow 
velocities  at  corresponding  elevations  are  much  lower  in  the  light-breaking 
case  (Fig.  19).  Such  a distinction  night  not  be  evident  from  Dean's 
stream  function  analysis,  since  it  is  derived  from  a surface  profile. 

Thus,  again,  the  measure  of  breaking  intensity  is  shown  to  be  the  (growth) 
rate  at  which  a wave  deforms  as  it  approaches  limiting  steepness. 


*The  droplets  shown  are  actually  higher  velocity  data  separating  from  the 
wave  face. 


PHASE  VELOCITY  284  cm/sec 
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Velocity  field  within  a lightly  breaking  0.66  hz  wave  (E-4,  Fig.  9)  at  the  same  pre 
breaking  stage  as  that  in  Fig.  16. 
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The  principal  conclusions  resulting  from  flow  velocity  measurements 
can  be  summarized  as  follows. 

o Owing  to  lateral  symmetry,  the  flow  fields  within  unsteady  progres- 
sive waves  of  steepness  lower  than  about  0.10  can  be  calculated  with 
reasonable  accuracy  by  Stokes'  5th-order  steady  wave  theory. 

o As  opposed  to  a singular  cuspated-crest  instability  (Price,  1971), 
breaking  appears  to  be  the  culmination  of  a continuous  process  whose 
effects  extend  well  below  trough  level.  Beginning  at  a steepness, 
H/L  = 0.1,  it  is  characterized  by  increasing  profile  asymmetry, 
tilting  of  velocity  isolines,  and  progressive  increase  of  velocities 
at  all  elevations,  until  the  forward  crest  face  erupts  as  a jet, 
slightly  beneath  maximum  crest  elevation,  and  within  which  flow 
velocities  attain-  or  slighly  exceed-  crest  speed.  Visible  e' idence 
of  crest  instability  occurs  well  in  advance  of  this  last  stage. 

o The  distinction  between  light  and  heavy  breaking  is  principally  con- 
fined to  the  relative  magnitudes  of  flow  velocities,  which  increase 
with  growth  rate. 


5.5  Surface  Profiles 


Of  all  wave  parameters,  surface  elevation  is  the  easiest  to  measure 
in  the  laboratory,  if  not  the  field.  All  later  consideration  of  potential 
energy  and  energy  loss  are  based  on  such  measurements.  In  this  section, 
the  steepness  limit  to  which  higher  order  theory  reasonably  describes  wave 
profiles  is  examined,  together  with  a more  sensitive  measure  of  asymmetry 
for  waves  of  higher  steepness. 


5.5.1  Low  Steepness  Profiles 

Figure  15  shows  that  the  velocity  isolines,  flow  directions,  and 
surface  profile  of  a young  wave  (B-3,  Fig.  9 and  Table  I)  all  exhibit 
lateral  symmetry.  Figure  20  compares  the  complete  profile  of  this  same 
wave  to  those  of  the  same  steepness,  ll/I,  = 0.10,  computed  from  Stokes' 
theory  to  3rd-order  (Kinsman,  1965,  p.  251),  5th-order  (Beach  Erosion  Board 
1941),  and  15th-order  (Monkmeycr  and  Kutzbach,  1965).  All  profiles  were 
normalized  to  L*,  and  the  measured  still-water  level  was  matched  to  the 
mean  of  the  computed  profiles.  Although  there  appears  to  be  a significant 
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improvement  in  going  from  3rd  to  5th  order,  the  distinction  between  5th 
and  15th  order  is  smaller  (1  cm)  than  the  measurement  uncertainty  (2  cm). 
Thus,  as  with  phase  speed  and  flow  velocity,  the  5th-order  solution  can 
be  regarded  as  adequate  to  describe  wave  profiles  for  H/L  s 0.10. 


- q-2  Asymmetry  in  Steep  Waves. 


As  shown  in  Figs.  16,  17,  and  19,  regardless  of  growth  rate,  profiles 
for  waves  steeper  than  0.10  become  increasingly  asymmetric.  These  changes 
are  most  pronounced  near  the  wave  crest  and,  in  heavily  breaking  waves, 
occur  very  rapidly.  Figure  21  shows  two  consecutive  photographic  profiles 
of  an  0.66  hz  wave  0.156  sec.  apart  (0.12  wavelengths)  as  it  approaches 
breaking.  Within  this  short  interval,  the  near-vertical  crest  face  has 
developed  into  a plunging  jet,  within  which  flow  velocities  exceed  the 
phase  velocity  by  about  8%. 

Because  growth  rate  is  difficult  to  assess  under  field  conditions,  a 
search  was  made  for  a quantitative  measure  of  asymmetry  which  might  be  cor- 
related with  growth  rate  as  an  alternative  parameter.  Of  several  possibil- 
ities, the  third  moment  of  slope  A(x,t)3  appeared  to  give  most  consistent 
results.  For  waves  moving  in  the  positive  x-direction,  this  moment  is 
given  by: 

x+L/2 

A(x,t)3  = - i'^L-  f (8n(x ' ,t) /9x ’ ) 3dx'  (4) 

x-L/2 

y 

where  -fl,"  is  an  arbitrary  normalization.  From  previously-demonstrated 

stat ionarity  (Section  5.2),  the  space-time  equivalent  of  (4)  is,  approximately, 

t +'!'/  2 

A(x,t)3  = - 1/T  / [9n(x,t  Vot'J-Mt'  (5) 

t-T/2 

If,  for  any  time  increment,  -T/2<t"<T/2,9n(t+t")/9t  = -9n(t -t”) / 3 1 , 
the  profile  is  symmetrical  and  A(x,t)  =0.  If  the  forward  wave  face  is 
steeper  than  the  back  face,  A(x,t)>0,  and  conversely.  In  computer  format, 

(5)  was  further  approximated  as: 

A(x,t)3  = 1/N  >-t  (nj-nj_i)3/At2  (6) 

where  nj  is  the  ith  elevation  sample  in  a wave  staff  time  scries  of  sample 
interval,  At. 
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Fig.  21 


HORIZONTAL  DISTANCE,  x (cm) 


Photographic  profiles  of  heavily  breaking  0.66  hz  wave  over  a time 
interval  of  0.156  sec.,  showing  rapid  acceleration  of  crest  to  form 
high  velocity  jet. 
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Table  III  lists  the  asymmetries  calculated  from  (6)  for  the  first 
three  or  four  uniform  waves  averaged  over  about  20  repetitive  runs  for 
light  breaking  series  E-4  and  E-5  (Fig.  9),  and  heavy  breaking  series  B-3, 

9,  and  10,  together  with  their  growth  rates  and  similarly-averaged 
steepnesses.  Although  there  appears  to  be  a consistent  increase  in  average 
asymmetry  between  young  and  pre-breaking  waves,  the  asymmetry  variations 
among  individual  waves  of  a series  are  substantially  larger  than  corresponding 
steepness  variations,  and  effectively  preclude  any  distinction  between 
heavy  and  light  breaking,  or  between  pre-breaking  and  breaking.  This 
(essentially  negative)  result  suggests  that  steepness  is  still  a more  re- 
liable index  of  breaking,  and  that  the  possibility  of  determining  growth 
rate  from  surficial  evidence  appears  doubtful. 
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5. 6 Wave  Energy 


One  of  the  principal  objectives  of  these  experiments  was  to  quantify 
energy  lost  by  wave  breaking,  and  to  parameterize  energy  loss  rates  in 
terms  of  frequency,  steepness,  and  growth  rate.  Because  of  decreasing 
measurement  accuracy  in  flow  velocities  below  about  50  cm/sec.,  quantitative 
estimates  of  kinetic  energy  were  restricted  to  wave  crests  above  trough 
level,  and  then  only  up  to  the  breaking  point.  ,‘onsequent  ly,  all  energy 
estimates  reported  here  refer  to  potential  energies  calculated  from  the 
variance  of  surface  elevation  averaged  over  time,  converted- where 
necessary-to  space  averages  by  virtue  of  demonstrated  stationarity  (Section 
5.2).  Energy  estimates  within  the  breaking  zone  suffered  from  corresponding 
inaccuracies  owing  to  the  difficulty  of  precisely  determining  the  instantaneous 
surface  elevation,  and  must  be  regarded  only  as  upper  bounds.  Despite 
these  shortcomings,  by  assuming  that  that  total  mean  energy  density  was 
equally  partitioned  between  potential  and  kinetic,  and  that  total  energy 
flux  was  conserved  before  breaking,  it  was  possible  to  construct  reasonable 
and  consistent  representations  of  the  energy  equilibria  prevailing  in  growing 
and  breaking  waves  (Section  5.7). 
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TABLE  III 

ASY^TRY  VS.  STEEPNESS 


f = 0.66  hz 


Station 

Order 
of  Wave 
Arrival 

Growth 

Rate 

G 

Steepness 

H/L 

Asymmetry 

A 

(cm/sec) 

B-3 

1 

0*250 

0.093 

0.04±0. 18 

2 

tt 

0.095 

0. 14±0. 19 

(young) 

3 

it 

0.097 

-0.25±0.06 

4 

it 

0.095 

-0.25±0.16 

E-4 

1 

0.266 

0.125 

0.77±0.03 

(pre- 

2 

tt 

0.123 

0.60±0.04 

breaking) 

3 

tt 

0.125 

0.32+0.27 

4 

tt 

0.132 

0.73±0.05 

E-5 

1 

0.295 

0.131 

1.24+0.18 

2 

tt 

0.138 

0.90±0.01 

(breaking) 

3 

tt 

0.129 

1.17+0.03 

4 

it 

0.132 

0.47±0.09 

B-9 

1 

0.535 

0.118 

0.77+0.16 

(pre- 

2 

tt 

0.134 

1.51±0.21 

breaking) 

3 

tt 

0.132 

0.66+ 

B-10 

1 

0.592 

0.132 

0.84+0.06 

2 

tt 

0.136 

1.33±0.09 

(breaking) 

3 

tt 

0.134 

0.61+0.04 
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5.6.1  Energy  Balance  Within  a Convergent  Channel 

In  the  most  general  sense,  the  energy  balance  equation  for  surface 
waves  on  deep  water  with  no  mean  flow*  is  given  by  (Phillips,  1969,  p.  50): 

3E(x,t)/3t4*^t(x,y,t)  = Q(x,t)  , (7) 

where  E(x,t)  is  the  total  energy  density,  ?t(x,y,t)  the  total  energy  flux, 
and  Q(x,t)  the  energy  dissipation  per  unit  area.  The  symbol  ^ is  the 
horizontal  gradient  operator,  i3/3x+j3/3y,  where  i and  j are  unit  vectors 
in  the  x and  y directions,  respectively.  In  this  study,  paddle  stroke  was 
kept  at  a constant  amplitude  after  an  initial  growth  period,  so  that 
3E/3t  = 0,  £t(x,y,t)  = ?t(x,y),  and  Q(x,t)  = Q(x) . 

Figure  22  shows  a section  of  the  tapered  channel  from  above.  Arrows 
denote  unit  vectors  n normal  to  the  sides  of  the  quadrilateral  area  A. 

Let  p be  the  perimeter  of  A.  If  equation  (7)  is  integrated  over  A and 
Green's  theorem  is  applied  to  the  integral,  one  obtains 

pf  $**t(x,y)  d\  = p/  £t(x,y)-n  dp  = Q(x)  dA.  (8) 

Letting  ?t(x,y)«n  * 0 along  the  channel  walls,  and  neglecting  the  y-component 
and  y-dependence  of  ^t(x,y)  and  Q(x)  along  the  open  ends  of  A introduces  an 
error  of  order  1 - cos  6°-20'  £ 0.01.  Then,  if  the  left  end  of  the  quadri- 
lateral is  a distance  Xj  from  the  paddle  and  the  right  end  a distance  x2, 

(8)  reduces  to 

X1 

w(x2)Ft(x2)  - w(x1)Ft(x1)  = / Q(x) w(x)  dx,  (9) 

where  Ft(x)  is  the  x-component  of  Ft(x,y). 

Equation  (9)  expresses  the  conservation  of  energy  flux  in  a slowly 
converging  channel.  It  can  be  put  in  more  useful  form  incorporating  the 
growth  rate,  G = -L*(3w/3x)/w,  by  holding  x^  constant,  and  differentiating 
with  respect  to  x-variable: 

-Ft(x)G/L*+3Ft(x)/3x-Qv(x)-Qr(x)  = Qb(x) , (10) 

where  Qy,  Qr>  and  are  losses  attributable  to  viscous  dissipation 
reflection,  and  breaking,  respectively.  The  problem  of  determining  break- 
ing losses,  then,  reduces  to  evaluating  the  terms  on  the  LHS  of  (10). 

*In  these  experiments,  data  sampling  was  concluded  before  transport  could 

establish  even  transient  mean  flow. 
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5.6.2  Energy  Flux 

For  small  amplitude  waves  propagating  in  one  direction,  energy  flux 
is  conventionally  given  by: 

Ft(x)  = CgE(x),  (11) 

where  C = 3w/3k  = g/4nf  is  the  group  velocity,  and  E(x)  the  total  energy 
S 

density.  Hayes  (1973)  has  recently  shown  that  (11)  also  applies  to  finite 

amplitude  waves  in  deep  water,  provided  that  the  proper  dispersion  relation 

is  used  in  determining  C , and  that  E(x)  is  equally  partitioned  between 

S 

potential  and  kinetic  energy. 

For  present  purposes,  group  velocity  was  determined  from  (2b)  by 
taking  f = 2-nr/cu  and  = u/k,  and  solving  for  the  differential  relation: 

Cg5  = 3w/3k  = (g/4rf) [1  + |{kH)2  + ~ (kH)4]  (12) 

Although,  as  previously  stated,  total  kinetic  energy  was  not  determined  in 
these  experiments,  it  is  shown  in  Appendix  G that  the  error  in  assuming 
energy  equipartition  in  the  steepest  and  most  asymmetric  waves  observed  in 
the  converging  channel  is  substantially  less  than  4%.  Accordingly, 

pgE(x)  = 2pgEp(x),  (13) 

where  pgEp(x)  is  the  potential  energy  density. 

Potential  energy  density  was  calculated  here  from  the  discrete  time 
series: 

t+T/2  7 

PgEn(x,t)  = £ n(x,f)^t’  (14) 

p t-T/2 

which  can  be  compared  with  its  spatial  equivalent  (Kinsman,  1965,  p.  146): 

x+L/2 

pgEn(x,t)  = Hf  Z n(x’,trAx'.  (15) 

P ^ x-L/2 

Their  normalized  difference,  [Epx(x,t)-Ep(x,t)]/Ep(x,t) , is  estimated  to 
be  less  than  0.05  (Appendix  C) . 

Combining  equations  (12),  (13),  and  (14),  the  expression  for  total 
energy  f’ux,  F‘t(x),  becomes, 

Ft(x)  = 2F(x)  = 2pgCg5Ep(x).  (16) 


(16) 


I 


c 


I 
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5.6.3  Viscous  Dissipation 


The  spatial  degradation  of  wave  amplitude  by  viscous  dissipation  within 
layers  at  the  bottom,  sides  and  free  water  surface  of  waves  in  a channel  is 
given  by  (Van  Dorn,  1965) : 


i2  = HjCXp-  / 


X2  2 , , . ^ l/2r,.kw(x) (l+cosh^kD)+sinh  2kD) 


w(x) 


(v/4irf)‘^[(i 


2kD+sinh  2kD 


)]k  dx,  (17) 


where  v is  the  kinematic  viscosity,  and  k(x)  = 2nf/Cr.  Energy  dissipation 

b 2 

was  calculated  from  this  equation  by  approximating  energy  density  by  n • 
This  was  an  accurate  estimate  because  the  spectral  peak  at  the  paddle 
frequency,  f,  contained  over  90  percent  of  the  wave  energy  (Appendix  B) . 

Over  the  travel  distances  listed  in  Table  II,  maximum  energy  losses, 
computed  before  breaking  for  all  frequencies,  varied  from  3 to  4%,  which 
was  small  compared  to  the  energy  increase  (about  100%)  occuring  during 
wave  growth.  Accordingly  it  was  assumed  that  Qv  = 0 in  Eq.  (10). 


5.6.4  Wave  Reflection 


Of  several  possible  approaches  to  estimating  the  amount  of  wave 
energy  reflected  from  the  converging  channel  walls,  that  adopted  here 
essentially  reduces  to  calculating  energy  flux  with  the  assumption  of  no 
reflection  losses,  and  then  estimating  the  attendant  error  from  the  variance 
of  the  computed  results  from  lines  of  conserved  flux.  Energy  flux  is 
discussed  in  Section  5.7.1,  and  the  above  calculations  are  described  in 
Appendix  E.  The  results  indicate  that  less  than  one  per  cent  of  incident 
energy  was  reflected  before  breaking,  so  that  Q = 0 in  Eq.  (10). 


5.6.5  Potential  Energy  After  Breaking 

Because  breaking  turb'ii f-"ce-particularly  on  the  forward  wave  face- 
could  be  expected  to  render  continuous  wave  staff  records  somewhat  ambiguous, 
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independent  evidence  was  sought  for  obtaining  more-accurate  estimates  of 

mean- square  elevation  after  breaking.  Figure  23  shows  photographically 

that  the  crests  and  troughs  of  even  a heavily  breaking  wave  were  relatively 

smooth,  so  that  wave  height  could  be  estimated  within  about  2 cm. 

Figure  24  is  a plot  of  mean-square  surface  elevations  before  and  after 

2 

breaking  versus  the  Stokes  linear  wave  energy  parameter,  H /16,  for  two 
frequencies  and  a wide  range  of  steepnesses.  The  solid  curve  for  Stokes 
5th-order  potential  energy  was  computed  from  the  expansions  given  by 
Mason  (Beach  Erosion  Board,  1941),  using  k = 2wf/Cg.  From  this  figure, 
it  can  bo  concluded  that: 

o The  breaking  and  non-breaking  distributions  are  reasonably  colinear. 
with  increasing  scatter  at  high  steepness. 

o The  fact  that  mean- square  elevations  among  breaking  waves  do  not 

differ  significantly  from  those  among  non-breaking  waves,  independent 
of  frequency  or  steepness,  has  possible  application  to  energy 
determinations  from  measurements  of  surface  elevation  in  storm  seas. 

o Both  data  distributions  fall  systematically  beneath  the  Stokes 
5th-order  curve,  by  amounts  that  increase  with  steepness. 

The  dashed  curve  in  Fig.  24  was  computed  from  the  empirical  expression: 

ij  2 i|2 

E = £I_(1-A!i-),  (18) 

P 16^  ,2 

where  A = 6.70  is  the  solution  of  (18)  for  H = H and  E = E = mean- 

* p p 

square  elevation  for  a limit  Stokes  wave. 

While  this  curve  may  have  no  physical  significance,  it  does  pass 
through  the  point  defined  by  the  potential  energy  for  a Stokes  wave  of 
limiting  steepness,  and  it  proviaos  a convenient  and  adequate  description 
of  the  data  distribution. 


5 . 7 Energy  Flux  and  Breaking  Losses 

Having  reasonably  established  the  energy  flux  equilibrium  within 
the  convergent  channel  in  terms  of  mcasureablc  variables,  it  remains  to 
attempt  to  relate  local  equilibria  to  frequency,  steepness  and  growth 
rate,  in  order  to  quantify  energy  loss  rate  as  a measure  of  breaking 
intensity. 


Fig.  23  Trough  (upper)  and  crest  (lower  plate)  of  a 0.66  iiz, 
heavily  breaking,  wave  18  m from  the  paddle.  The 
vertical  scales  are  in  decimeters. 
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and  reflected  energy  can  be  neglected,  so  that  Qv  = Qr  = = 0,  energy  flux 

is  conserved,  and  equation  (10)  can  be  written: 

G/L*  = -(3w/3x)/w  = OFtx/3x)/Ftx  (19) 

Integrating  (19)  between  the  limits  xq  and  x,  and  recalling  (18),  one 
obtains: 

U/Go  « Fx(x)/Fx(x0),  (20) 

which  shows  that  energy  flux  is  proportional  to  growth  rate  when  normalized 
to  common  initial  conditions.  For  present  purposes,  Fq  = Fx(xq)  was 
arbitrarily  chosen  to  be  that  flux  corresponding  to  G = Gq  = 0.19,  since, 
for  all  wave  series  in  Table  II,  there  happened  to  be  a station  where  growth 
rate  was  close  to  this  value,  and  no  wave  had  yet  broken. 

Evidence  that  the  experimental  non-breaking  energy  flux  determinations 
are  consistent  with  (20)  is  shown  in  Figs.  25-27  for  all  wave  series  at  the 
three  operating  frequencies  listed  in  Table  II.  Since  the  data  series  were 
ordered  in  terms  of  increasing  paddle  amplitude,  commencing  with  strokes 
so  small  that  no  breaking  occured  within  the  test  section  (Series  A) , 
these  data  never  depart  significantly  from  the  conserved-flux  lines, 
whose  slopes  are:  F/FqGo  = 1/0.19  = 0.67.  With  increasing  wave  amplitude 

(initial  st 'epness) , breaking  of  successive  wave  series  - as  indicated  by 
dashed  branchings  from  these  lines  - occured  with  decreasing  intensity  at 
correspondingly  lower  growth  rates.  If  this  qualitative  presumption 
is  true,  the  branch  points  marking  the  onset  of  breaking  energy  loss 
should  also  be  correlated  with  growth  rate.  The  difficulty  here  is  the 
precise  determination  of  the  branch  points.  For  a given  wave  scries, 
normalized  energy  flux  is  essentially  the  product  of  the  square  of  wave 
steepness  (Fig. 24)  multiplied  by  group  velocity,  both  of  which  vary  mono- 
tonically  with  growth  rate,  so  that  their  product  has  no  well-defined 
maximum  at  breaking.  Accordingly,  the  best  estimates  of  breaking  inception 
were  here  taken  as  the  dashed  intercepts  of  the  post-breaking  flux  lines 
in  Figs.  25-27  with  the  lines  of  conserved  energy  flux.  These,  normalized 

3 

b>  L*f  to  remove  frequency  dependence,  are  plotted  against  growth  rate  in 


GROWTH  RATE,-L  ^/w 

*dx 


Fig.  25  Energy  flux  vs.  growth  rate  at  f = 0.66  hz  for  seven  initial  wave 
steepnesses. 
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Fig.  26  Energy  flux  vs.  growth  rate  at  f = 0.73  hz  for 
five  initial  wave  steepnesses. 
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Fig.  27  Energy  flux  v 
wave  steennes 
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Fig.  28.  Expectedly,  the  data  distribution  much  resembles  that  for  breaking 
steepness  (Fig.  10),  exhibiting  considerable  scatter  with,  perhaps,  a 
slight  tendency  to  increase  with  growth  rate.  The  dashed  line,  drawn  through 
the  data  by  inspection,  is  convex  upward  in  accordance  with  the  intuitively- 
reasonable  assumption  that  normalized  energy  flux  must  approach  some  asymptotic 
limit  as  G + “. 


5 . 7. 2 Breaking  Energy  Loss  Rate 


Subsequent  to  breaking,  Qb  i 0,  and  (10)  becomes: 


3Ft/3x-FtG/L*  = Qb.  (21) 

2 

Recalling  (16),  and  making  use  of  the  operator,  3/3x  = G (3/3G)/L*,  (21) 
can  be  written: 

pgQb  - 2P^f-  fF0  ^ (ergs/cm2-sec)  , (22) 

which  must  be  satisfied  along  the  dashed  branched  lines  in  Figs.  25-27. 

A more  convenient  measure  of  energy  loss  rate  is  its  ratio,  per  wave  cycle, 
to  the  average  energy  present,  E = 2Fx/C^,.: 


% 

fE 


= _G  fa 


G 3Fy  1 
F„  SG^l 


(23) 


In  this  equation,  all  quantities  are  known,  except  for  the  normalized 
flux  gradient  parameter,  (1/F0) (3FX/3G) , which  is  given  by  the  slopes 
of  the  above  branch  lines. 

As  with  the  branch  points,  the  quantitative  determination  of  potential 
flux  gradients  is  somewhat  subjective,  since  they  depend  upon  small 
variations  in  group  velocity  and  potential  energy  amid  significant  data 
scatter.  Qualitatively,  the  dashed  linos  in  Figs.  25-27  appear  to 
have  systematically  negative  slopes  at  low  transition  growth  rates  (large 
initial  wave  amplitude  and  weak  breaking);  slope:,  tend  to  increase  with 
transition  growth  rate  (greater  breaking  intensity),  and,  perhaps,  even 
go  slightly  positive  during  the  early  stages  of  intense  breaking. 

Despite  these  initial  (transitional)  tendencies,  all  slopes  tend  toward 
zero  if  breaking  proceeds  long  enough. 
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Fig.  28  Energy  flux  at  breaking  inception  vs.  growth  rate. 
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As  regards  equation  (23),  the  physical  interpretation  to  be  placed 
on  these  tendencies  is  tentatively  the  following.  A zero  flux  gradient 
implies  that  breaking  energy  loss  is  always  in  equilibrium  with  the  energy 
input  imposed  by  channel  convergence,  the  rate  of  both  being  determined 
by  a wave's  initial  steepness  and  growth  rate.  Just  prior  to  breaking 
transition,  there  appears  to  be  tendency  for  all  waves  to  overshoot  their 
respective  equilibrium  breaking  flux  rates,  and  to  decay  back  to  the 
equilibrium  rate  as  breaking  continues.  Both  overshoot  and  decay  rate 
appeal  to  increase  with  breaking  intensity. 
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For  lack  of  sufficient  data  to  properly  quantify  these  transient 
approaches  to  equilibrium,  breaking  energy  losses  were  instead  estimated 
by  considering  the  extreme  ranges  of  relevant  variables: 

o potential  energy  flux  gradient;  -2<(1/Fq) (3Fx/3G)<0 

o breaking  steepness;  0. 125<H/L<0. 140 

o normalized  group  velocity;  0.64<C^/C*<0.70 

Figure  29  is  a plot  of  normalized  energy  loss  rate  at  breaking  transition 
(Eq.  23)  versus  growth  rate.  The  straight  lines  bounding  the  "error"  bars 
give,  roughly,  the  uncertainty  attached  to  determinations  of  energy  flux 
rate  for  all  wave  series  in  Table  II.  Uncertainty  of  group  velocity  is 
given  by  the  shaded  bands  centered  on  these  lines.  The  upper  line  might 
be  considered  the  best  estimate  of  breaking  overshoot,  and  the  lower 
line  that  for  equilibrium  breaking  loss,  corresponding  to  Q^/fn  = -0.93G 
and  -0.67G,  respectively.  Thus,  within  the  range  of  these  experiments, 
breaking  loss  rates  amount,  roughly,  to: 


I 


% 

Breaker 

Growth  Rate 

Q,  /fl:  (per  cent/cycle) 

ft 

Type 

G = -L*(3w/3x)/w 

Transition 

F.qui  librium 

1 

Light 

0.2 

19 

13 

£ 

% 

z 

Heavy 

0.6 

57 

39* 

*This  value  is  an  extrapolation,  since  no  equilibrium  breaking  was  observed 
at  this  growth  rate  because  of  limited  channel  length. 

The  above  rates  indicate  substantial  energy  losses--evcn  for  light  breaking, 
and  the  median  value  (26%)  for  equilibrium  breaking  is  about  3000  times 
larger  than  that  per  mean  wave  cycle  estimated  for  sea  waves  (Longuct- 


Uiggins,  1969b).  While  this  result  may  have  little  significance  to  natural 
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conditions,  the  concept  of  growth  rate  is  perfectly  general.  Its  application 
to  the  more  natural  case  of  biharmonic  waves  breaking  as  a result  of  con- 
structive interference  is  discussed  in  the  next  section. 


6.  DISCUSSION 

This  section  is  devoted  first  to  a brief  summary  and  interpretation 
of  the  principle  results  of  this  study  and,  second,  to  an  extension  of 
these  results  to  an  estimate  of  the  energy  loss  rate  attributable  to 
waves  breaking  solely  by  virtue  of  constructive  interference  among  components 
of  a narrow  spectrum. 


6 . 1 Results  Summary  and  Interpretation 


The  growth  towards  breaking  of  uniform  trains  of  periodic  deep  water 
waves  propagating  within  a straight-sided  convergent  channel  was  found  to 
be  characterized  by  three  successive,  steepness-dependent  phases: 

1)  a young  phase  (0/<H/L<0. 10) , wherein  surface  profiles  and  internal 
flow  fields  are  symmetric  and,  together  with  phase  (crest)  speeds, 
are  roughly  in  accord  with  those  piedicted  by  Stokes  steady  wave 
theory  to  5th  order; 

2)  a pre-breaking  phase  (0. 10<H/L<0. 125) , wherein  there  is  no 
obvious  surface  manifestation  of  crest  instability,  but  where 
crest  profiles  are  markedly  unsymmetric,  flow  isolincs  are  con- 
cave upwards  and  tilted  in  the  direction  of  wave  propagation, 
and  strong  velocity  gradients  develop  just  beneath  the  crest; 

3)  a break ing  phase  (0. 1 25<H/L<0. 14) , involving  verticality  of  the 
forward  crest  face,  obvious  crest  instability,  and  increased 
tilting  and  curvature  of  sub-crest  fluid  velocity  iso’incs  so  as 
to  form  a plunging  jet,  within  which  flow  velocities  equal  or 
exceed  the  local  crest  speed. 

Except  for  phase  speed,  which  continues  to  accord  with  Stokes  5th- 
order  theory,  none  of  the  properties  described  in  phases  (2)  and  (5),  above, 
are  predicted  by  any  steady  wave  theory.  The  fact  that  their  degree  and 
speed  of  occurence  increases  with  growth  rate,  G,  suggests  that  they  in- 
volve unsteady  (accelerative)  processes.  At  the  same  time,  evidence  for 
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continued  stationarity  (space-time  equivalence)  right  up  to  breaking 
suggests  that  these  changes  are  not  specific  to  a convergent  channel,  but 
would  occur  similarly  in  a parallel  channel,  if  its  width  could  somehow  be 
mechanically  decreased  at  corresponding  rates. 


An  alternative  and,  perhaps,  intuitively  clearer  picture  cf  the 
internal  flow  field  changes  associated  with  terminal  wave  growth  is  given 
by  the  "streamline"  representations  in  Fig.  30a,  b,  c,  whose  indices 
correspond,  respectively  to  the  three  growth  phases  described  above,  as 
well  as  to  the  flow-isoline  representations  for  a heavily-breaking  0.66hz  j 

wave  shown  in  Figs.  15-17.  In  a coordinate  system  moving  at  wave  speed,  j 

streamlines  can  be  defined  as  curves  everywhere  parallel  to  the  local  j 

flow  direction  whose  spacing  varies  inversely  with  flow  velocity.  While  f 

the  curves  shown  were  subjectively  contoured  from  observational  flow  and  j 

profile  data1,  they  provide  a qualitative  insight  to  breaking  kinematics.  J 


In  a steady  inviscid  (young)  wave  (Fig.  30a),  the  surface  is  a stream- 
line, and  all  subsurface  streamlines  are  homologous  nested  curves,  symmetric 
about  the  wave  crest.  With  increasing  steepness,  (Fig.  30b) , profile 
asymmetry  is  relected  by  corresponding  streamline  asymmetry,  at  least  to 
trough  depth.  However,  with  breaking  (Fig.  30c),  the  flow  is  no  longer 
steady;  streamlines  can  terminate  in  the  surface,  or  inviscid  flow  can 
create  a vortex.  While  the  vortex  shown  is  speculative,  the  indicated  flow  j 

reversal  corresponds  to  fluid  velocity  exceeding  phase  speed,  and  offers  j 

a possible  explanation  for  the  observation  that  the  breaking  jet  always  j 

issues  somewhat  below  the  wave  crest.  If  this  concept  has  physical  j 

validity,  one  expects  that  energy  dissipation  could  arise  both  from  vortex  j 

shear  and  jet  advection  out  of  the  wave. 

As  regards  the  energy  balance  in  growing  waves,  the  following  points  ! 

can  be  made.  The  fact  that  potential  energy  flux  was  found  to  be  conserved  j 

right  up  to  breaking  transition  is  perhaps  the  best  indication  that 
potential  flow  persists,  and  that  kinetic  energy  flux  is  similarly  con- 
served. This  conclusion  is  supported  oy  the  arguments  for  energy  equi- 
partition  presented  in  Appendix  G.  The  onset  of  breaking-defined  by  a 
relatively  abrupt  drop  in  the  local  energy  flux  rate-was  invariantly 


*An  attempt  was  made  to  calculate  streamlines  directly  by  least-squares 
fitting  the  observed  flow  fields  to  7th-ordcr  polynomials.  But  even  this 
order  was  insufficient  to  follow  details  of  crest  flow. 


associated  with  a transient  normalized  energy  loss  rate  of  about, 

Q^f/E  = -0.93G,  that  decayed  within  a wavelength  to  an  equilibrium  level 
of  about  two-thirds  this  rate.  Such  an  equilibrium  implies  that  waves  lose 
energy  through  breaking  at  just  the  rate  at  which  it  is  supplied  by  channel 
convergence.  Thus  the  equilibrium  energy  flux  is  also  proportional  to 
growth  rate.  Since  the  former  is  approximately  the  product  of  group 
velocity  times  twice  the  local  energy  density,  both  of  which  are  weak, 
monotonical ly  increasing  functions  of  steepness  (Eqs.  12  and  18),  the 


above  proportionality  implies  that  equilibrium  breaking  steepness  must 
also  increase  weakly  with  increasing  growth  rate.  Although  the  experi- 
mental data  are  not  sufficiently  accurate  to  quantify  this  premise,  the 
entire  range  of  breaking  steepness  (0. 125<H/L<0. 14)  indicates  that  a 
steepness  increase  of  only  12  percent  triples  the  breaking  loss  rate. 
(13-39  per  cent/cycle). 


| 
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6.2  Interactive  Breaking 

A logical  extension  of  these  results  is  the  question  of  their 
application  to  the  breaking  of  wind  waves  at  sea.  Despite  the  recognized 
randomness  of  sea  waves,  when  the  wind  is  steady-as  in  the  trades-the 
surface  is  often  characterized  by  the  appearance  of  rather  regular  wave 
groups.  Owing  to  small  directional  differences  of  energy  propagation, 
the  groups  are  staggered  in  a 'chicken-wire'  pattern.  Under  such  circum- 
stances, Donelan,  et  al  (1972)  have  observed-and  confirmed  by  simple 
experiment-that  breaking  occurs  systematically  at  the  group  antinodes  at 
intervals  of  about  twice  the  period  of  the  principal  waves.  While 
breaking  of  this  type  requires  a fairly  narrow  frequency  and  direction 
spectrum,  there  is  increasing  evidence  (Rye,  1974)  that  the  storm  wave 
spectrum  is  rather  narrower  than  that  earlier  proposed  by  Pierson  and 
Moskowitz  (1964).  Lastly,  preliminary  experiments  with  interactive 
breaking  among  waves  of  a unidirectional  biharmonic  scries  in  the  present 
(non-convergent)  wave  channel  fully  substantiates  the  premise  that 
regular  breaking  can  be  produced  by  constructive  interference  in  the 
i.jsence  of  wind.  The  ensuing  discussion  advances  the  hypothesis  that 
variable  steepness  and  growth  rate  may  account,  in  part,  for  the 
qualitative  features  observed. 
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Consider,  as  the  simplest  example,  a linear  biharmonic  wave  train  \ 
composed  of  two  sinusoids  having  the  same  amplitude,  a,  but  different  i 
radian  frequencies,  nj>n9,  and  wave  numbers  k^>k7,  related  by  k = n /g:  ! 

n(x,t)  = a cosCkjX-njti  + a cos(k?x-n2t) . (24)  l 

Although  breaking  may  depend  weakly  upon  nonlinear  effects,  the  time  and  \ 
place-and  perhaps  the  duration-of  breaking  should  coincide  roughly  with  ■ 
this  linear  model.  ■ 


The  general  nature  of  breaking  interaction  is  better  visualised  by 
an  alternative  form  of  (24): 


n(x,t)  = 2a  cos(k  x~n  t)cos(k  x-n  t) , (25) 

where  the  subscripts  c and  m represent  mean  sums  and  differences,  respective- 
ly, of  the  component  frequencies  and  wavenumbers  defined  above.  Observed 
in  space  at  a fixed  time  (Figure  31a)  or  in  time  at  a fixed  point,  the 
wave  form  is  a carrier  wave  of  amplitude  2a.  frequency  nc,  and  wavenumber 
k , modulated  bv  a sinusoidal  envelope  of  frequency  n , and  wavenumber  k . 

In  a tank  experiment  where  such  a disturbance  is  generated  by  a wave  paddle, 
the  carrier  waves  are  those  actually  seen  by  an  observer;  they  move  at 
"phase  velocity",  C£  - nn/kc,  passing  through  the  invariant  modulation 
envelope,  which  propagates  a "group  velocity",  C’m  = nj:^.  At  any  instant, 
the  number  of  carrier  waves  in  a group  is,  N = kc/2k  . 

By  analogy  to  the  convergent  channel,  it  is  now  supposed  that  carrier 

wave  growth  is  everywhere  specified  by  frequency,  n , and  the  local 

values  of  steepness  (11^,/L^)  and  growth  rate,  Gc,  and  that  breaking  will 

occur  whenever  the  local  steepness  exceeds  some  critical  value,  say, 

II  /L  >0.13. 
c c 

The  parametric  transformation  is  most  simply  developed  in  a 

coordinate  system  moving  at  phase  velocity,  Cc>  so  that  x = x’  + C^.t, 

where  x’  is  the  space  coordinate  fixed  to  some  particular  wave.  If  the 

time  coordinate  is  redefined  as  t = t'-fk  x'/n  B)  , where  B = (k  n -k  n )/k  n 

me  m c t m c e 

and  t'  = ('  at  the  antinode  of  the  modulation  envelope  (Fig.  31a), 
equation  (25)  can  be  rewritten  as: 

n(x',t')  = 2a  cos(k  x ' )cos(  Bn^.t  ' ) . (26) 

Thus,  as  a carrier  wave  passes  through  a node  of  the  modulation  envelope, 
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Fig.  31  (a)  Biharmonic  wave  packet  composed  of  carrier  waves 

modulated  by  a low  frequency  envelope, 
fb)  Equilibrium  breaking  of  duration  t^  reduces  post- 
breaking envelope,  as  shown  by  dashed  curves. 


-78- 


it  grows  with  time,  t',  and  its  local  steepness  is  given  by: 


H /L  = 11  k /2tt  = ~ cos  (Bn  t')  = ~[1  + tan2(Bn  t’)]~1/2  (27) 

CCCC  L c l c 

c c 

The  carrier  wave  energy  flux,  F , is  proportional  to  its  potential 
energy  density,  hence  to  the  square  of  the  modulation  amplitude: 

F 'vcos2(Bn  t*)  . (28) 

c c 

Accordingly,  by  analogy  to  equation  (19),  local  growth  rate,  C:c>  can  be 
defined  as  the  normalized  time  rate  of  change  of  Fc: 

Gc  = ~ * jr  f£r  = tan(Bncf)  (29) 

c c 


Evidently,  Gc  decreases  from  infinity  at  the  modulation  nodes  (where  the 
carrier  amplitude  vanishes)  to  zero  at  the  antinodes,  and  has  intermediate 
values  within  the  range  (0.6>G>0.2)  of  the  convergent  channel  experiments. 


That  value  of  growth  rate,  G^,  corresponding  to  the  onset  of  breaking 
(Hc/Lc  * 0.13)  can  be  determined  through  substitution  of  (27)  into  (29): 


G*. 

c 


— 4 71  B [ ( 


2a/Lc 

0.065 


)2-l] 


1/2 


(30) 


which,  expectedly,  is  a function  only  of  the  carrier  steepness  index, 

2a/L,  and  the  frequency  ratio  parameter,  B.  If,  as  in  the  convergent 
channel,  the  subsequent  breaking  history  is  characterized  by  equilibrium 
between  growth  rate  and  breaking  energy  loss  rate,  breaking  steepness  will 
thereafter  remain  constant  until  the  wave  reaches  the  antinode  of  the 
original  modulation  envelope,  (G  = 0),  where  breaking  ceases  as  growth 
rate  begins  to  diminish.  Thereafter,  the  reduced  modulation  envelope  is 
defined  by,  r>  = 0.131.  cos(Bn  t * ) , as  shown  bv  the  dashed  curves  in  Fig.  31b. 
Similarly,  the  breaking  duration,  tj  , in  carrier  wave  periods,  Tc,  can 
be  determined  from  (29)  and  (30) : 


, ....  1 _ r,2a/l,c,2  .,1/2 

tb/  c " ‘ 2uB  drctjn  0 . 065  " J 


(31) 


The  mean  fractional  breaking  energy  dissipation  rate  per  carrier 

wave  period  is  here  most  conveniently  determined  by  considering  the  wave 

group  as  a black  box  whose  input  energy  flux  is  C Ej  and  output  flux  is 

C E0.  The  fractional  dissipation  rate,  O , is  then  the  ratio  of  the  energy 
m i l 


-79- 


lost  in  a wave  period  to  the  average  energy  present: 


% ‘ WE,  - E2)/LcEl  • (VCc»>  - 


(32) 


Since,  after  breaking,  the  energy  density  ratio  is  F^/Ej  = [0. 065/(2a/Lc) ] ' 
and  recalling  thatCm  and  C^are  functions  only  of  the  biharmonic  frequency 


ratio  a = n^/n^,  (32)  becomes,  finally: 


% = 


1 + a 
(1  + a) 


0.065,2 
2a/L  ' 


(33) 


All  of  the  foregoing  results  are  summarized  in  Eig.  32,  in  which 
isolines  a, a,  G^.,  and  Q are  drawn  in  a coordinate  system  with  steepness 
index  as  ordinate  and  breaking  duration  as  abscissa.  While  it  would  have 
been  more  orthodox  to  use  a as  an  independent  variable  than  t^/T^  the 
above  representation  provides  clearer  distinction  between  the  other  variables. 
The  reasonable  range  of  natural  sea  conditions  is  effectively  bounded  by 
the  lines  0.6<a<0.9,  corresponding,  roughly,  to  N = 1 and  N = 5 waves  per 
wave  group,  as  indicated  by  numbers  in  parentheses.  The  shaded  region 
covers  the  range  of  experimental  values  of  these  parameters  determined 
for  the  convergent  wave  channel. 


The  conclusions  inferred  from  this  (admittedly-simplified)  model  for 
biharmonic  breaking  are,  tentatively,  the  following.  Note,  first,  that 
breaking  losses  are  determined  almost  solely  by  the  steepness  index,  and 
are  almost  independent  of  breaking  kinematics.  This  is  a result  of  the 
assumption  of  equilibrium  breaking  within  the  black  box,  which  uniformly 
truncates  the  wave  envelope,  figure  32  indicates  that  for  small  a (wide 
component  frequency  separation),  there  will  be  only  one  or  two  waves  in  a 
group;  growth  rate  will  be  high  and  breaking,  presumcably,  intense,  but 
breaking  duration  is  short.  Conversely,  for  large  a,  there  will  be  many 
waves  in  a group,  associated  with  low  growth  rate  and  weak  breaking  of 
longer  duration.  Within  the  region  of  practical  interest,  these  effects 
tend  to  compensate,  so  that  total  breaking  loss  remains  nearly  constant. 


Second,  as  regards  application  of  the  biharmonic  model  to  wind  wave 
breaking,  honguct -Higgins  (1969bJ  and  Hasselmann  (1974)  indirectly  arrive 
at  estimates  of  order  10  ^ and  10  , respectively  for  breaking  energy 


i A 1 1 breaking  loss  curves  in  Eig.  32  drop  hyperbolical ly  to  zero  near  the 
left  margin.  This  is  because  breaking  is  limited  by  decreasing  steepness 
when  there  is  less  than  one  wave  per  group. 


STEEPNESS  INDEX, 2a/Lc 


BREAKING  DURATION,  tb/Tc 

I so  lines  of  growth  rate  Cd  at  breaking,  frequency  ratio  a = n?/n 
and  equilibrium  breaking  energy  loss  Qc>  as  functions  of  earner 
wave  steepness  index  2a/ 1.  and  breaking  duration  tj,/Tc.  Shaded 
area  shows  region  of  overlap  between  present  experiments  and 
narrow-band  ocean  wave  i ntcrac' ions. 
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loss  per  mean  wave  cycle  within  the  equilibrium  range  of  the  energy 
spectrum.  The  first  estimate  supposes  a narrow  spectrum,  and  the  second 

a breaking  duration  short  compared  with  the  'period'  of  the  principal  waves. 

-2 

With  these  restrictions,  Fig.  32  indicates  loss  rates  of  order  10  '. 

However,  the  present  model  is  unidirectional  and  assumes  infinite  crest 
length.  The  geometric  extension  to  bidirectionality  would  have  the  effect 
of  additional  crest-wise  modulation  of  wave  steepness,  and  hence  breaking. 

A rough  estimate  suggests  that  this  would  reduce  loss  rates  by  about  an 
order  of  magnitude. 

Thus,  despite  neglect  of  randomness,  the  influence  of  wind  shear, 
and  internal  energy  transfer  among  frequencies,  the  present  biharmonic 
model,  extended  to  include  directionality,  appears  to  be  not  too  far 
removed  from  reality  to  warrant  further  study.  The  principal  results 
suggest  that  breaking  losses  depend  mainly  upon  steepness  and  breaking 
duration  (via  Gp  and  are  largely  independent  of  frequency.  The  former 
variables  are  amenable  to  field  observation  and  laboratory  study,  and  form 
the  basis  for  further  work  in  this  important  and  relatively  unexplored 
area  of  wave  dynamics. 
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APPENDIX  A 

Determination  of  Fluid  Velocity  and  Direction 

The  functional  part  of  the  D1SA  55A81  hot  film  probes  used  in  these 
experiments  consists  of  a wedge-shaped  quartz  rod  about  3 mm  in  diameter 
on  whose  chisel  edge  is  sputtered  a resistive  nickel  film  element  1 mm 
long  by  0.35  mm  wide.  When  connected  to  a DISA  55D01  constant  temperature 
anemometer,  DISA  claims  a frequency  response  of  over  30  khz  in  water  over 
the  velocity  range  10<U<2500  cm/scc.  This  response  is  sufficiently  rapid 
to  assure  equilibrium  output  by  the  time  the  probe  tip--initially  in  air-- 
has  penetrated  a distance  of  0.1  mm  into  a vertical  wave  front  moving  at 
300  cm/sec. 

The  anemometer  output  is  an  analog  voltage  V,  functionally  related 
to  the  relative  flow  velocity  U past  the  probe  tin  by: 


V2  - v2  = A(0)Un 


(A- 1) 


where  Vq  is  the  no-flow  voltage,  0 is  the  probe  >aw  angle  (Fig.  A- la)  . 

The  physical  "constant"  n and  the  normalized  directional  coefficient,  A(0) , 

2 

are  determined  for  each  probe  by  calibration.  ' Figure  A-2  shows  a 
representative  set  of  calibration  curves  made  in  the  500  ft  Lockheed 
towing  tank  in  San  Diego,  and  covering  the  range  of  flow  velocities 
(3CKU<550  cm/sec)  and  yaw  angles  f 0< 0<tt / 2 ) common  to  these  experiments.* 

The  following  general  comments  can  be  made. 

o The  fact  that  these  curves  are  not  straight  lines  shows  that  n r 
constant,  but  decreases  slightly  with  increasing  velocity.  Because 
fitting  these  curves  with  even  the  simplest  polynomial  approximation 
greatly  increases  the  complexity  of  determining  angular  response, 
they  were,  instead,  approximated  by  two  straight  lines  of  slope  n 
for  30<U<65  cm/sec,  and  n,  for  <>5<U<300  cm/sec,  respectively. 

As  later  shown,  this  simplification  introduced  errors  of  less  than 
6 cm/sec  over  the  entire  velocity  range. 


^Sincc  an  investigation  of  pitch  response  indicated  less  than  1 of 
error  due  to  cross  channel  velocities,  the  pre-breaking  measurements 
were  essentially  two-dimensional. 

^Whilc  the  angular  response  of  hot  wi re  probes  can  be  approximated  from 
physical  principles  Oiinzc,  1959,  (*.h.  2),  flow  ’orturbat  ions  induced  by 
a film  probe's  shape  and  attitude  proscribe  si.  , 1c  analysis. 

•'The  calibrations  actually  extended  over  yaw  angles  from  -r/2  to  ~/2 
in  order  to  resolve  small  probe  nssymmet rics. 


I ; o The  fact  that  all  curves  have  every  where  the  same  slope  at  any 
I given  velocity  made  it  possible  to  use  the  same  values  of  na  and 
I . n^  for  all  yaw  angles. 

I o Repeat  calibrations  at  different  temperatures  in  a closed  water 
I tunnel  at  the  University's  Applied  Mathematical  and  Engineering 
| Science  Laboratory  established  the  temperature  independence  of 


na,  and  n^,  and  yaw  response. 

Since  polar  plots  of  normalized  yaw  response  g(9)  = (A(6)  - A(n/2)]/ 
[A(0)  - A(it/2)]  exhibited  slight  sinuous  perturbations,  the  response  of 
each  probe  was  expanded  in  truncated  Fourier  series: 

6 

g(0)  = Z g cos  2i0.  (A-2) 

i = 0 

The  Fourier  coefficients  g^  were  determined  for  13  yaw  angles  in  ^-in- 
crements and  for  7-10  velocities  within  the  range  50<U<550  cm/sec.  by 
standard  methods.  Representative  values  of  for  a single  probe, 
averaged  over  seven  velocities,  together  with  their  respective  standard 
deviations  Ag.,  are  listed  in  Table  A-l.  For  the  first  two  coefficients, 
which  account  for  about  90  percent  of  the  total  contribution  to  g(6) , 
the  proportional  error  Ag^g^  is  less  than  0.01.  Such  small  deviations 
indicate  that  the  probe's  angular  response  is  essentially  independent  of 


velocity. 

Table  A-l 

index  (i) 

0 

1 

2 3 

4 

5 

6 

gi 

0.649 

0.442 

-0.122  0.043 

-0.020 

0.016 

-0.008 

&8i 

. 005 

.004 

.005  .003 

.003 

.005 

.002 

The  directional  ambiguity  A0(0)  associated  with  these  deviations  Agj 

can  be  estimated  by  setting  ^ ^ Ag . + A9  = 0 and  solving  for 

i = 0 Ri  1 30 

A9(0) : 


A9(0)  = Ag.cos2i9/  £ 2ig.sin2i0 

i = 0 1 i = 0 1 


(A- 3) 


which  can  be  evaluated  from  the  data  in  Table  A-!  for  any  yaw  angle  9. 

The  results  are  shown  in  Table  A-2,  from  which  it  is  apparent  that  the 
directional  ambiguity  of  a single  probe  is  always  small,  except  near  . 

9=0,  for  which  sin  2 i 0 -•  0 and  the  angular  response  becomes  indeterminate. 
However,  for  a pair  of  orthogonal  probes,  the  directional  ambiguity  is 
always  determinable  and  small  (see  below). 


.1  mi*. 
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Table  A- 2 


Flow  angle  0 (degrees) 
Ambiguity  0 (degrees) 


10  20 
4.0  0.3 


30 

0.6 


40 

0.3 


50 

0 


60 

0.1 


70 

0.1 


80 

0.1 


90 

2.0 


Having  verified  the  form  of  the  probe  response  equation  (A-l) , it 
remains  to  generalize  the  factor  A( 0)  in  terms  of  the  normalized  directional 
response  g(0)  and  some  measureable  response  parameter,  say  A(tt/4) , since 
all  probes  were  dynamically  calibrated  at  0 = tt/4  in  these  experiments.  To 
do  this,  it  was  found  convenient  to  define  an  additional  directional 
constant  which,  like  the  coefficients  g^  could  be  evaluated  from  a complete 
directional  calibration  of  each  probe,  and  stored  for  later  use.  Of 
several  alternatives,  the  definition:  g*  = gQ  + A(tr/2)/[(A(0)  - A(tt/2)], 

where  gQ  is  the  first  Fourier  coefficient  in  A-2,  leads  to  a particularly 
simple  expression  for  directional  response: 


Ate)  =A(»/2)[f^i-!^] 


(A-4) 


The  values  of  g*  obtained  for  each  probe  were  all  found  to  fall  within 
the  range  1.0<g*<1.15,  and  can  be  thought  of  as  being  related  to  slight 
physical  eccentricities.  The  simplification  g*  = 1.0  would  lead  to 
directional  errors  of  a few  degrees,  which  is  sufficient  for  most  practical 
purposes.  However  the  constants  actually  determined  were  employed  in  all 
data  analysis. 

Since  the  anemometer  output  is  the  product  of  a velocity  function 
and  a directional  function,  two  probes  were  needed  to  resolve  these 
functions  separately.  In  the  present  experiments,  the  two  probes  were 
mounted  in  the  vertical  plane  of  the  probe  turret  so  that  their  tips 
nearly  abutted  at  right  angles  (Fig.  A-lb).  Although  the  turret  was 
rotatable  in  45°  increments,  consider  first  the  configuration  shown,  for 
which  the  unknown  flow  has  a velocity  U and  a direction  making  an  angle 
0 with  the  (horizontal)  bisector  of  the  probe  axes.  If  the  probe 
(anemometer)  outputs  are  defined  by:  P^  = (V2  - V2)  j = UniAj(0j)  = UniAj(0  + tt/4) 

and  P2  = (V2  - V2)2  = Un2A2(02)  = Un2A2(0  - tt/4),  then,  for  all  0 within 
the  range  tt/4>0>-tt/4,  the  ratio  of  the  directional  coefficients  is  unique. 

If  the  probes  had  identical  velocity  responses,  so  that  n^  = n2»  0 could 
be  obtained  by  computing  Pj/P,,  anc*  searching  a table  of  Aj/A2  for  the 
corresponding  value.  However,  even  where  n^  / n2,  one  can  use  an  analogous 


‘j^smims^&miM 


procedure.  Letting  m = t^/n^,  we  can  write: 

Pr/P2  = Aim(0  + */4)/A 2(0  ' */4)>  (A_5) 

from  which  0 can  be  found  from  a table  of  A^/A^  Spurious  ratios  of 
Pim/P2,  such  as  may  occur  with  one  probe  out  of  water,  will  fall  outside 
the  range  of  this  table.  Having  found  0,  the  corresponding  flow  velocity 
U is  given  explicitly  for  either  probe  by: 


U = exp  [ C 1/n) ln(P/A) ] 


(A-6) 


Separate  directional  calibration  of  individual  probes  in  the  Lock- 
heed towing  tank  was  considered  to  be  much  more  efficient  and  inherently 
more  accurate  than  calibrating  them  by  pairs  in  the  turret  head,  which 
would  have  required  a much  larger  number  of  test  runs.  Once  all  probes 
were  calibrated  and  their  Fourier  coefficients  g.  and  directional 
constants  stored  on  punch  cards,  they  could  be  ex- .banged  in  the  turret 
at  any  time  when  there  was  evidence  of  malfunction.  Altogether,  only 
five  different  probes  were  used  in  the  course  of  three  months  of  data 
taking.  Probe  failure  was  usually  evidenced  by  progressive  loss  of 
sensitivity  after  several  hundred  hours  of  use. 

During  data  runs,  the  turret -mounted  probe  pairs  were  dynamically 

calibrated  at  the  fixed  yaw  angles  0^  = -Sj  = tt/ 4(0  = 0)  over  the  velocity 

range  80<U<550  cm/sec,  as  follows.  The  probe  turret  mounting  strut  was 

first  elevated  to  an  indexed  position  and  rotated  90°  until  the  turret 

locked  into  place  with  the  probe  bisector  vertically  centered  on  the  0.9  cm- 

diameter  discharge  orifice  in  the  bottom  of  the  jet  calibrate  cylinder 

(Fig.  A- 3).  This  cylinder  was  about  10  cm  in  diameter  and  2 m high,  and 

was  equipped  with  an  orifice  plug  valve  (initially  closed),  and  11  pairs 

of  shorting  contacts  up  its  side.  The  cylinder  was  filled  with  water 

pumped  directly  from  the  wave  channel  so  as  to  provide  a calibration  at 

the  in  situ  temperature  of  the  channel.  The  plug  valve  was  then  opened, 

and  the  issuing  jet  impinged  symmetrically  upon  the  probe  tips  about  1.2  cm 

beneath.  As  the  water  level  in  the  cylinder  fell  consecutively  below 

each  contact  pair,  an  electrical  signal  instructed  the  computer  to  sample 

two  anemometer  voltages  Vj  and  The  turret  strut  was  then  returned  to 

the  vertical  and  lowered  until  the  turret  was  immersed,  and  the  computer 

instructed  by  push  button  to  sample  the  no- flow  anemometer  outputs 

V , and  V The  computer,  having  previously  been  supplied  with  jet 
ol  o2 


Fig.  A-3  Method  of  mounting  strut  for  hot-film  velocity  probes  for 
elevation,  reversal,  and  calibration. 
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1 2 2 
velocities  corresponding  to  the  11  contact  pairs  , then  plotted  ln(V  - V ) 

vs  In  U for  each  probe.  These  points  were  fitted  with  a straight  line 

and  its  intercepts  at  In  U = 4 and  6 read  off,  from  which  the  calibration 

constants  n and  A(ir/4)  were  computed  and  stored  for  later  use.  The  data 

points  plotted  in  Fig.  A-2  show  how  closely  a typical  cylinder  calibration 

fits  the  45°  yaw  angle  curve  obtained  from  the  directional  towing  tank 

tests.  From  analysis  of  randomly  selected  cylinder  calibrations  it  was 

concluded  that  the  computed  directional  resolution  was  accurate  to  within 

2°  at  most,  and  usually  much  less  than  1°,  and  that  the  velocity  resolution 

above  65  cm/sec  was  within  3%.  The  constant  n for  the  velocity 

range  30<U<65  cm/sec  could  not  be  determined  from  cylinder  calibrations. 

However,  repeated  spinning  tank  tests  at  0 = 0 showed  that  all  probes 

had  very  similar  low-velocity  responses,  and  rms  value  n = 0.54  was  used 

to  compute  all  velocities  within  this  rai.ge.  The  best  error  estimate  for 

low  velocities  was  3°  in  direction  and  46  cm/sec  in  velocity.  This 

amounted  to  fairly  large  percentage  errors  in  velocity,  but  because  the 

velocities  were  low,  the  errors  did  not  significantly  affect  the  overall 

results  or  the  conclusions  drawn. 


s 


\ 


^These  were  determined  independently  by  calculation  and  by  weighing  the 
discharge  at  constant  contact  levels. 
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maker  produces  particle  displacements  closely  akin  to  those  of  the  desired 
waves  in  stationary  sequence,  free  secondary  waves  will  be  generated  in' 
phase  with  the  second  harmonic  of  the  paddle  frequency.  Because  the  primary 
and  secondary  waves  travel  at  different  speeds,  the  wave  profile  may  suffer 
progressive  distortion.  This  phenomenon  is  particularly  important  in 
shallow  water  (Van  Dorn,  1975),  but  was  examined  here  as  a possible  i 

contributor  to  wave  asymmetry  even  though  the  depth  was  always  half  the 
primary  wavelength. 

I 

Figure  B-l  shows  three  frequency  (potential  energy)  spectra 
calculated  by  Fourier  methods  (Jenkins  and  Watts,  1968)  from  25-second 

i 

wave  staff  records  of  0.66  hz  waves  at  stations  D-2,  F-2,  and  G-2  j 

(Table  D) . Although  relatively  close  to  the  wave  paddle  (593  cm),  the 

waves  were  sufficiently  high  (h  = 50  cm)  and  steep  (H/L  = 0.12)  as  to  ! 

represent  extreme  conditions  for  all  wave  series  in  these  experiments. 

All  three  spectra  are  coherent  to  at  least  the  4th  harmonic,  with 
sharp  peaks  centered  on  integer  multiples  of  0.66  hz.  The  primary 

i 

peak  contained  at  least  70%  of  the  potential  energy  in  a band  0.04  hz  j 

wide,  and  over  90%  in  a 0.2  hz  band.  The  2d  harmonic  contains  only 

about  6%  of  the  energy  of  the  primary.  Corresponding  spectra  at  in- 
creasing distance  from  the  paddle  exhibited  ro  detectible  phase  shift  of 
the  2d  or  higher  harmonics.  Accordingly,  it  was  concluded  that  free 
waves--if  present — were  too  low  to  significantly  affect  the  results 
reported  here. 


from  a time  series  of  surface  elevation  at  a fixed  location  (Eq.  14)  is 

equivalent  to  that  which  would  be  obtained  from  a spatial  profile  at  a J 

corresponding  instant  (Eq.  15) . An  estimate  of  the  associated  error  was  % 

JE 

obtained  as  follows.  -5| 

X 

i 

Surface  elevation  can  be  expressed  as  a sum  of  harmonics:  | 

n(t,x)  = aj(x)cos0  + a2(x)cos20  + a3(x)cos30  + ..,  (C-l)  -J 


where  k = 30/3x,  u = 30/3t,  x is  distance  along  the  channel  from  an 
arbitrary  origin,  and  k is  a slowly  varying  wavenumber.  If  variation  is 
small,  an(x),  n = 1,..,  can  be  represented  by  the  first  two  terms  of  a 
Taylor's  series. 


an(x)  — + a^x  + ^ = 2,...  (C-2) 

The  integral  definitions  of  equations  (14),  and  (15)  are  respectively. 


1 


T/2 


— / 
2T  -T/2 


1 

2L 


L/2 


-L/2 


n(t,0) 2dt, 
2 

n(0,x)  dx. 


Using  equations  (C-l),  (C-2),  and  the  identities 
2 2 3 2 

and  y cos  ny  dy  = ti  /3  + it/2n  , where  n = 1,2,3,. 
written: 


(C-3a) 


(C-3b) 

* 2 

/ ycos  ny  dy  = 0, 
these  can  be 


F = — E a2 
p 8tt  n0’ 

n 

_ 1 r r 2 .22  f7T  1 ,, 

px  8 it  an0  anl  6 4Ttn2^  ’ 

n 

and  their  normalized  difference  is 


e„,-£p  L'l°-6  * °-54  »2i 


I 


10 


20 


(C-4a) 


(C-4b) 


(C-5) 


I 

I 


Since  most  of  the  wave  energy  is  contained  in  harmonic  peaks  of  the 


This  table  is  a computerized  printout  of  wave  heights  and  mean- 
square  surface  elevations  calculated  automatically  from  time  series 
elevation  measurements  at  individual  wave  staffs.  The  table  is  divided 
into  three  sections  in  order  of  increasing  paddle  frequency,  f.  Each 
section  is  subdivided  into  consecutive  wave  series,  A,B,C...,  in  order 
of  increasing  paddle  stroke.  Each  series  listing  gives,  from  left  to 
right,  a station  number,  paddle  distance,  channel  width,  and  the  corres- 
ponding wave  heights  and  potential  energies  of  the  first  four  uniform 
waves  averaged  over  three  consecutive  data  runs  at  a wave  staff  9 cm  from 
the  channel  sidewall.  Corresponding  data  for  an  additional  cross-channel 
staff  at  some  stations  appear  in  the  last  two  columns. 

Wave  heights  were  calculated  as  the  maximum  elevation  differences  in 
the  time  series,  and  energies  as  half  the  time-variance  of  elevation, 
averaged  over  a wave  period.  These  values  were  arithmetically  averaged 
over  the  triplicated  repeat  runs  to  give  the  table  entries. 
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APPENDIX  E 


Wave  Reflection 


Systematic  deviations  of  F/Ffl  from  the  line  of  conserved  energy 
flux  can  be  used  to  calculate  an  upper  limit  on  the  reflection  term, 
Qr(x),  in  equation  (10).  Wave  elevation  time  series,  n(x,t),  can  be 
represented  as  a sum  of  reflected,  R,  and  transmitted,  I,  Fourier 
coefficients: 


I|(x.t)  = £ ivy  * 

i 1 1 

where  wavenumber  k = 27rf/C_  (p.  27).  Since  the  phase  angle,  0., 

3 J 

determines  the  relative  complex  phases  of  reflected  and  transmitted 

coefficients,  both  1^  and  R^are  real.  Total  potential  energy  is 


(E-l) 


-Lgn2  = |pg  [I"?  + R?  + 2I.R.cos(k.x  - 0.)], 


(E-2) 


and  by  equation  (16) , the  total  energy  flux  is  2CgSpgn  . Since  over 
90  percent  of  this  sum  is  contained  in  the  frequency  band  centered 
on  the  paddle  stroke  frequency,  f,  the  energy  flux  is  approximated  by  that 
single  term.  The  normalized  potential  energy  flux  is  approximately 

EI2/Iq]  + [2IRcos(kx  - 6}/Iq], 

where  L and  R are  the  transmitted  and  reflected  Fourier  coefficients  at 
the  paddle  frequency,  respectively,  and  IQ  is  a transmitted  Fourier  co- 
efficient corresponding  to  the  normalizing  energy  flux,  F_.  It  it 

2 2 2 2 U 
assumed  that  Rp«I0,  and  R «I  . An  upper  limit  on  normalized  energy  flux 

fluctuations  due  to  the  last  term  in  (E-2),  above,  can  be  estimated  from 

the  10  percent  maximum  deviations  from  the  lines  of  conserved  energy 

flux  observed  in  Figures  25-27,  at  stations  before  breaking.  Station 

spacing  was  close  enough  so  that,  for  at  least  half  uf  the  stations, 

2 

|cos(kx-0) | >0.5,  which  leads  to  0. l>IR/In>R/I . Therefore  the  ratio  of 

2 2 u 

reflected  to  transmitted  energy,  R /I  was  less  than  0.01. 
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APPENDIX  F 
Error  Estimates 

Error  estimates  associated  with  velocity  measurements  are  discussed 
in  Appendix  A.  This  section  deals  with  the  maximum  probable  errors  in 
the  operational  data  of  waveheight,  steepness,  potential  energy,  and  energy 
flux  listed  in  Tables  I and  II  (Section  S.l),  which  were  derived  by 
averaging  from  Table  D,  which,  itself,  represents  averages  from  multiple 
data  runs. 

Except  for  steepness  and  energy  flux,  which  involve  phase  and  group 
velocities,  repectively,  all  errors  in  the  above  parameters  arise  from 
variances  in  the  measurement  of  surface  elevation,  compounded  by  averaging 
between  cross-channel  wave  staffs,  and  among  four  uniform  waves  in  a train 
and  from  three  to  15  repetitive  data  runs. 

While  the  number  of  individual  staffs,  uniform  waves,  and  repeat 
runs  was  always  too  small  to  apply  usual  tests  of  statistical  normality, 
the  smallness  of  their  fractional  standard  deviations  (-0.04)  implies 
that  their  meaned  distributions  are  considerably  more  likely  to  be  normal 
(Hoel,  1947,  p.  69).  Accordingly,  with  this  assumption,  if 
(Table  D)  are  ensemble  averages  of  wave  height  over  N repeat  runs, 
including  cross-channel  staffs,  let  H = H./2N  and  H = H/4  be  the  cross- 


.rr  2,1/2 


channel  and  four-wave  means,  respectively,  and  a(H)  = [<H  > - <H>  ] 
and  o(H)  = a(H)/2  be  their  respective  standard  deviations  (Topping,  1960). 
If  similar  terminology  applies  to  the  mean-square  elevations  (potential 
energies,  E^)  in  Table  D,  then  their  respective  95%  confidence  limits 
are,  approximately,  H ± o(H)  and  ± o(Ep) . 

Defining,  now,  the  corresponding  fractional  errors  in  wave  height  and 
potential  energy  as,  a(H)/H  and  o(Ep)/Ep,  using  the  previously  estimated 
fractional  error  for  phase  velocity,  oC_/Cr  = 0.03  (Sect.  5.1.1),  and 
assuming  that  for  group  velocity  C^g  to  be  no  larger,  the  respective 
error  estimates  for  steepness,  H/L  = H/CgT,  and  for  energy  flux, 

Ft  = 2CggEp,  become,  finally  (Ibid.,  p.  82): 


,77.2 


= me  + 0,032]1/2; 


(P-1) 


H/L 
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and 


sffitl  _-  [SgnL  ♦ 0T032]1/2. 
Ft  Ep 


(P-2) 


The  range  of  these  estimates  varied  slightly  between  stations  and 
wave  series,  but  the  overall  95%  confidence  limits  for  fractional  errors 
listed  in  Section  5.1  are  considered  to  be  conservative. 
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APPENDIX  G 


Energy  Partition  in  an  Inhomogeneous  Wave  Field 


For  any  conservative  wave  system,  total  energy  is  the  sum  of 
potential  and  kinetic  energy;  for  small -amplitude  waves,  the  two  are 
equal.  However,  kinetic  energy  could  not  be  determined  in  the  present 
experiments  because  flow  velocities  were  measured  only  under  the  wave 
crests.  Since  no  means  of  calculating  kinetic  energy  with  inhomogeneous 
waves  of  finite  amplitude,  have  been  reported,  resort  was  made  to  a result 
obtained  by  Davies  (1951),  for  symmetric  waves  of  finite  amplitude.  To 
fifth  order,  he  gives: 


1.080 


(G-l) 


where,  A = 0.992  for  a Stokes  wave  of  limiting  steepness.  Subject,  then, 
only  to  the  uncertainty  imposed  by  inhomogeneity,  the  present  assumption 
that  Ep  = Ek  is  only  in  error  by  about  4%,  which  is  less  than  the  fractional 


error  estimate  (±  3%)  given  in  the  previous  section  for  the  cal- 


culation of  E itself. 
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